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Abstract 

A theory of transformation is presented for the diagonalization of a Hamiltonian that is quadratic in creation and annihilation 
operators or in coordinates and momenta. It is the systemization and theorization of Dirac and Bogoliubov-Valatin transfor- 
mations, and thus provides us an operational procedure to answer, in a direct manner, the questions as to whether a quadratic 
Hamiltonian is diagonalizable, whether the diagonalization is unique, and how the transformation can be constructed if the 
diagonalization exists. The underlying idea is to consider the dynamic matrix. Each quadratic Hamiltonian has a dynamic 
matrix of its own. The eigenvalue problem of the dynamic matrix determines the diagonalizability of the quadratic Hamilto- 
nian completely. In brief, the theory ascribes the diagonalization of a quadratic Hamiltonian to the eigenvalue problem of its 
dynamic matrix, which is familiar to all of us. That makes it much easy to use. Applications to various physical systems are 
discussed, with especial emphasis on the quantum fields, such as Klein-Gordon field, phonon field, etc.. 

PACS numbers: 02.10.UD, 02.30.Hq, 03.65.-w 

Keywords: Linear algebra. Ordinary differential equations, Quantum mechanics 



Contents 

I. Introduction 

A. Bogoliubov-Valatin transformation 

B. Equation of motion 

II. Diagonalization Theory of Bose Systems 

A. Dynamic matrix 

B. Necessary condition for diagonalization 

C. Sufficient condition for diagonalization 

D. Uniqueness 

III. Diagonalization Theory of Fermi Systems 

A. Existence 

B. Uniqueness 

IV. Application to Bose Systems 

A. The normal Hamiltonian 

B. The pairing Hamiltonian 

V. Application to Fermi Systems 

A. The normal Hamiltonian 

B. The pairing Hamiltonian 

VI. Generalized Bogoliubov-Valatin 
Transformation 

A. Theory of Dirac diagonalization 

B. Landau quantization 

C. Partial diagonalization 

VII. Field Quanta 

A. Klein-Gordon field 

B. Phonon field 

C. Dirac field 

VIII. Mathematical Essence of 
Diagonalizability 



6 

1 
J 



'Electronic address: ' xmwQnetra. nj u. edu. on | 



A. Heisenberg operator [42 

B. Time-polarized photons [43 

IX. Conclusions |45 
Acknowledgments |45 
References |45 

I. INTRODUCTION 

In the preface to the first edition of "The Principles of 
Quantum Mechanics" [l| , Dirac said: "The growth of the 
use of transformation theory, as applied first to relativity 
and later to the quantum theory, is the essence of the 
new method in theoretical physics." 

In this review, we shall present a theory of transforma- 
tion, which can perform diagonalization to the Hamilto- 
nian that is quadratic in creation and annihilation op- 
erators or in coordinates and momenta. Such transfor- 
mations appear most frequently in classical and quan- 
tum mechanics, statistical mechanics, condensed-matter 
physics, nuclear physics, and quantum field theory. 

For the sake of simplicity, let us begin with the so- 
called Bogoliubov-Valatin transformation. 

A. Bogoliubov-Valatin transformation 

In 1947, Bogoliubov Q introduced a novel linear trans- 
formation to diagonalize the quantum quadratic Hamil- 
tonian present in superfluidity. This method was later 
extended by Bogoliubov himself d i, i] and also by 
Valatin [1, 0] to the Fermi case in the theory of supercon- 
ductivity. It has ever since got widely used in different 
fields [a H J and known as Bogoliubov-Valatin (BV) 
transformation, including both the bosonic and fermionic 
versions. 



To show the underlying idea of the method due to 
Bogohubov and Valatin, let us consider the quadratic 
Hamiltonian, 

n ^ ^ 

H=Y1 {a^jc\cj + -7yclc] + 2^j»c,Cj), (1.1) 

where n > 1 is a natural number, and Ci and c\ are, 
respectively, the annihilation and creation operators for 
bosons or fermions. They satisfy the standard commuta- 
tion or anticommution relations. 



fermionic. Besides, the matrices a and 7 will not van- 
ish simultaneously; otherwise, the Hamiltonian H is zero 
trivially. 

If we define a new product between the two operators 
Ci (or c\) and Cj (or cj) as 



Cj, • C7 



i I J J 



(1.11) 



then Eqs. (|1.2p - (|1.4p can be expressed compactly as 



■0 • i/jt = /^^ 



where 
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= 5iJ^ 


(1.2) 




I 




Ci, Cj\ 


— Cj dz Cj Ci 


= 0, 


(1.3) 




±1 


1 




= clc] ± c]4 


= 0, 


(1.4) 


with / being the identity matrix of size n 



where 5ij is the Kronecker delta function. The coeffi- 
and 7y € C have the following symme- 



cients aij € 
tries, 



It, = T7j«, 



(1.5) 



where z* denotes the complex conjugate of z. Through- 
out this review, the complex field C will be used as the 
base field of the Hamiltonian H . 

Using the form of matrix, Eq. (|1.1[) can be written as 



i/-iv^MV'±itr(a), 



(1.6) 



where tr(A) denotes the trace of the matrix A. The V' is 
a column vector and i/j^ its Hermitian conjugate. 



V^t^ [ct, c] 



(1.7) 



where c and c) are the subvectors of size n. 



Cl 
C2 



'^^ [ C2, 



(1.8) 



Here A denotes the transpose of the matrix A. The co- 
efficient matrix M has the form, 



M ^ 



a 7 



(1.9) 



where a and 7 are the submatrices with and 7^^ as 
their their entries, respectively. Obviously, 



= a, 7 = =F7, — M. 



(1.10) 



That is to say, a and M are both Hermitian matri- 
ces whereas 7 is a symmetric or antisymmetric matrix, 
which is determined by whether the system is bosonic or 



(1.12) 



(1.13) 



To diagonalize the Hamiltonian of Eq. (|1.6p . Bogoli- 
ubov and Valatin introduced a linear transformation. 



Ad + BdK 



(1.14) 



where A and B are two square matrices of size n, and d 
and d^ are the vectors as follows, 



d = 



di 
d2 



d^^[di 4 



2' 



di 



(1.15) 



Here di and are the new annihilation and creation 
operators respectively, they satisfy the standard commu- 
tation or anticommution relations as in Eqs. (|1.2p - (|1.4p . 

which means. 



if ■ if 



where 



d 



[d\ d 



(1.16) 



(1.17) 



From Eqs. (fL7| . (fTTT)) and (fLM)) . it follows that 

^ = (1.18) 

where 



A B 
B* A* 



(1.19) 



Here A* denotes the complex conjugate of the matrix A. 
By the way, we note that such a form of T originates 
from the requirement that c and must be Hermitian 
conjugates of each other. For convenience, we shall call 
the operator vector such as -0 and Lp the field operator. 

Under the transformation of Eq. (jl.lSp . the Hamilto- 
nian of Eq. (|1.6p becomes 

H = ^(fi^'T^'MTip ± itr(a), (1.20) 
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where TtMT is the new coefRcient matrix. Meanwhile, 
Eq. (frT2|) turns into 



(1.21) 



where Eq. (|1.16p has been used. Obviously, this is a 
condition for the transformation of Eq. (jl.lSp . 

For the Hamiltonian H to be diagonalized with respect 
to the new annihilation and creation operators, it is nec- 
essary that the new coefficient matrix T^MT is diagonal, 
i.e., 



UJl 



LJ2 











W2„ 



(1.22) 



where uji for i — 1,2,-- - ,2n are the diagonal entries, 
they are real: lOi G R. Equation (|1.22p means that all 
the off-diagonal entries of the matrix T^MT must vanish 
identically. Under this condition, we have 



i=l i=l 



'^w„+,;±itr(a). (1.23) 



This is the so-called diagonalized form for the Hamilto- 
nian H . 

To sum up, Eqs. (|1.2ip and (|1.22p are the two condi- 
tions that must be fulfilled by the transformation matrix 
T. The former ensures the statistics of the system, i.e., 
the system will remains bosonic or fermionic after the 
transformation if it is bosonic or fermionic before the 
transformation, that is a physical requirement. The lat- 
ter ensures the diagonalization of the Hamiltonian, it is 
just a mathematical requirement. According to Bogoli- 
ubov and Valatin, the transformation matrix T can be 
determined from Eqs. (|1.21[) and (|1.22p . After the deter- 
mination of T , the diagonal entries oji for i = 1, 2, • • • , 2n 
will be obtained, which accomplishes the diagonalization 
procedure. That is the main idea of the Bogoliubov- 
Valatin transformation. 

As indicated by Eq. (|1.10p , the matrix M is Hermitian. 
So it can always be diagonalized by a unitary transfor- 
mation. At first glance, it seems as if the Hamiltonian 
of Eq. (|1.20p could be brought into diagonalization by 
the same unitary transformation. However, a close ob- 
servation shows that such a unitary transformation can, 
in general, neither take the form of Eq. (|1.19p nor meet 
the requirements of Eq. p.2ip although it always satis- 
fies the condition of Eq. (|1.22p . Therefore, the unitary 
transformation for the diagonalization of the coefficient 
matrix M can not generally diagonalize the Hamiltonian 
of Eq. (|1.20p . That is because both the field and 
the field are now the vectors of operators (quantum 
numbers) rather than the usual simple vectors of com- 
plex variables (classical numbers). For the latter, it is 
well known that a Hermitian quadratic form can always 
be diagonalized by the unitary transformation for the di- 
agonalization of its coefficient matrix. In short, the BV 



diagonalization for a quantum quadratic Hamiltonian is 
much more complicated than the unitary diagonalization 
for the usual Hermitian quadratic form of complex vari- 
ables. 

Finally, let us analyze the BV method in more detail. 
It can easily be seen from Eq. (|1.19p that the trans- 
formation matrix T has 4n^ independent unknown en- 
tries. However, Eqs. (|1.2ip and (ll.22p contain An? and 
An? — 2n constraints on T, respectively. That is to say, 
the constraints are much more than the total number of 
the free unknown entries of T. Therefore, there are two 
possibilities: (1) Those constraints are consistent with 
the requirement of T, and thus T has solutions. (2) The 
constraints are inconsistent with the requirement of T, 
and T has no solution. Theoretically, it is very difficult 
to judge which case will happen because, as indicated by 
Eqs. (|1.2ip and (|1.22p . the constraints constitute 8n'^ — 2n 
coupled quadratic equations for An? free unknowns. Fur- 
thermore, it will still be hard to solve for the multiple 
unknowns from the multiple equations of second degree 
even if there exist solutions for the matrix T. Mathemati- 
cally, these difficulties arise from the well-known fact that 
there is no much knowledge about the multiple equations 
of second degree with multiple unknowns at present. In 
practice, one often has to rely on experience and tricks 
when he uses the BV method to resolve practical prob- 
lems. 

To overcome those difficulties, we intend to develop 
a new theory for BV transformation. We expect that 
this theory can not only judge straightforwardly whether 
a quantum quadratic Hamiltonian is BV diagonalizable 
but also yield the required transformation by a simple 
procedure if the Hamiltonian is BV diagonalizable. That 
is the main objective of this review. 



B. Equation of motion 

As shown in the preceding subsection, the diagonaliza- 
tion scheme adopted by Bogoliubov and Valatin is merely 
algebraic. That is to say, the scheme treats the diagonal- 
ization just as a pure algebraic problem, it does not con- 
sider the physics in diagonalization at all. We would like 
to complement it with physical contents so as to find the 
necessary and sufficient conditions for the diagonalization 
of a quantum quadratic Hamiltonian. Simply speaking, 
we shall take into account the equation of motion of the 
system, i.e., the Heisenberg equation. 

To show the idea, let us consider the classical system of 
harmonic oscillators — the coun terp art of the Bose system 
with a quadratic Hamiltonian |11| . 

_^ 71 1 ^ 



i,3 = l 



(1.24) 

where qi and pi {i — 1,2,--- ,n) are, respectively, the 
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generalized coordinates and momenta, with q and p being 
the corresponding column vectors, 



(1.25) 



The K and V are the kinetic and potential matrices with 
Kij and Vij as their entries, respectively. They are both 
real and symmetric. 









pi 








P2 


q = 




, p = 






_ qn _ 




Pn 



K 
V 



K>0, 
V>0. 



(1.26) 
(1.27) 



It is worthy to emphasize that K is a positive definite 
matrix, that is because the kinetic energy is always pos- 
itive definite. In addition, the matrix V is only posi- 
tive semidefinite, the bottom of potential being chosen 
as zero. 

As well known, qi and (i = 1, 2, • • • , n) satisfy the 
following canonical relations. 



or equivalently. 



{Pi,Pj} 
{q^,Pj} 



q ■ q 

P-P 
q-p 



0, 
0, 

/. 



(1.28) 
(1.29) 
(1.30) 



(1.31) 
(1.32) 
(1.33) 



where {a, 6} denotes the Poisson bracket of a and b, and 
a ■ b = {a, b}. 

Of course, the Bogoliubov-Valatin scheme can 
be transplanted directly to diagonalize the classical 
quadratic Hamiltonian of Eq. (|1.24p with respect to the 
new generalized coordinates and momenta. However, we 
would rather here turn to another way — the canonical 
equation of motion. 

The canonical equation of motion can be deduced from 
the Hamiltonian of Eq. p.24p and the Poisson brackets 
of Eqs. pr^ - p:^ as follows. 



d 
di 
d_ 
dt 



q = {q, H} = Kp, 
p = {p,H} = -Vq. 



where t denotes the time. As a result, we have 



d2 

dt 



■^q = -KVq. 



(1.34) 
(1.35) 

(1.36) 



That is a homogeneous system of linear ordinary differ- 
ential equations with constant coefficients. 



From the theory of ordinary differential equations [12l | , 
we know that the solution of the homogeneous linear sys- 
tem above depends on the eigenvalue problem. 



LU^q = KVq. 



(1.37) 



This eigenvalue problem can be solved rigorously with 
the help of the Cholesky decomposition of K, 



K = 



(1.38) 



where Q is an invertible matrix. The existence of such a 
deconiposition stems mathematically from the positivity 
of K [lj|. By introducing a temporal variable ^, 



Eq. (|1.37p can be transformed into 



where 



= AC, 



A = QVQ = A > 0. 



(1.39) 
(1.40) 
(1.41) 



Just as V ^ the matrix A is still real, symmetric, and non- 
negative definite. So it can be orthogonally diagonalized. 



where 



SKS 



SS = SS= L 



T = 











(1.42) 
(1.43) 

(1.44) 



Here oj^ > {i — 1,2, ■ ■ ■ n) are the eigenvalues of A, and 
S the orthogonal matrix with the eigenvectors of A as its 
column vectors. 

From Eqs. ^jT^ . dUT]), (HH]), and it follows 

that 



where 



If we put 



T-^KVT = r. 



T = QS. 



T=[vt, V2, 



(1.45) 
(1.46) 
(1.47) 



where Vi (i — 1, 2, • • • ,n) denote the column vectors of 
T. Equation (|1.45|) shows that Vi are the eigenvectors of 
the matrix KV , 



(1.48) 
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belonging to the eigenvalues wf, respectively. In other 
words, they are the solutions of the eigenvalue problem 
of Eq. (|1.37p . Evidently, they are orthonormal and com- 
plete, 



where 



TGT 
TTG 



(1.49) 
(1.50) 



where G = K^^ . Namely, they constitute a n- 
dimensional Hilbert space with G as its metric tensor. 

Just as usual, the general solution of Eq. (|1.36p can be 
expanded in this Hilbert space as 



j{t) =Y,i^^{t>^ 



(1.51) 



where ipi{t) (i — 1,2,-- - ,n) are the expanding coeffi- 
cients. But, not as usual, we do not care here how to 
determine those coefficients from the initial conditions. 
Instead, we would rather view this expansion as a linear 
transformation, 



q{t) = Tilsit), 
where ^j{t) is the column vector. 



(1.52) 



Mt) 



(1.53) 



As will be seen later, this view is crucial for the diago- 
nalization of the Hamiltonian. Since T has full rank, the 
transformation is invertible. The inverse is 



(1.54) 



Physically, ■0(0 represents the new generalized coordi- 
nates, and q{t) the old ones. 

The corresponding transformation for the generalized 
momenta can be deduced from Eq. (|1.33p . As is well 
known, a Poisson bracket is a bilinear function of its two 
arguments. This together with Eq. (|1.33p indicates that 
there exists a duality relationship between p(t) and q{t) 
[T^ . This duality implies that p{t) will transform con- 
travariantly with q(t), i.e.. 



7rit)=Tp{t), 



(1.55) 



where 7r(i) represents the new generalized momenta. 

Under the transformation of Eqs. (|1.54p and (|1.55p . 
the Hamiltonian of the system becomes as follows. 



H 



1. 



2 ^ 

i=l 



1 



(1.56) 



^■ip = 0, 

TT • 7? = 0, 
■0-7? — I. 



(1.57) 
(1.58) 
(1.59) 



They are identical to the system of Eq. (|1.24p and Eqs. 
(|1.3ip - (|1.33p . with the Hamiltonian being diagonalized 
with respect to the new generalized coordinates and mo- 
menta. 

Proposition 1 A classical quadratic Hamiltonian such 
as Eq. {1.24-^ can be diagonalized with respect to the gen- 
eralized coordinates and momenta. 

This instance demonstrates clearly that the equation 
of motion is a very effective and powerful weapon for the 
diagonalization of a quadratic Hamiltonian, in compari- 
son with the method of the preceding subsection. We see 
that the equation of motion can generate a linear trans- 
formation in a very natural way, which can then not only 
diagonalize the quadratic Hamiltonian but also ensure 
the invariance of Poisson brackets. 

In the picture of diagonalization, the system is repre- 
sented by the normal modes of motion. To diagonalize a 
quadratic Hamiltonian is therefore equivalent to seeking 
the normal modes of the system. Of course, the natu- 
ral tool for seeking the normal modes is the equation of 
motion, from the point of view of physics. That is the 
physical interpretation for the diagonalization. All in all, 
the canonical equation of motion is a candidate way to 
the diagonalization of a classical quadratic Hamiltonian, 
in addition to the purely algebraic method due to Bogoli- 
ubov and Valatin. 

Since Heisenberg equation is the quantum counterpart 
of the canonical equation of motion in the classical me- 
chanics, it encourages us to try to employ Heisenberg 
equation to realize the diagonalization of the quantum 
quadratic Hamiltonian. That is the main idea of this 
review. 

Complying with this idea, we shall first study the BV 
transformation and diagonalization of the Bose system, 
we find that a complete theory can be established using 
the Heisenberg equation (Sec. [H]). And then we study 
the BV transformation and diagonalization of the Fermi 
system (Sec. IIIII) . it is parallel to the Bose case. Their 
applications are discussed in the following two sections 
(Sec. HV] and Sec. IV]). Afterwards, we would turn to 
studying the Dirac transformation and diagonalization, 
which concern coordinates and momenta. It is found 
that they are the generalizations of the BV transforma- 
tion and diagonalization (Sec. IVip . An advantage of the 
Dirac transformation and diagonalization is that they can 
be transplanted readily to the complex collective coordi- 
nates and momenta, and therefore have wide applications 
in field quantization (Sec. IVIip . Finally, we would like 
to clarify the mathematical essence of the transforma- 
tion and diagonalization. We find that the equation of 
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motion can be sublated. The transformation and diag- 
onalization have nothing to do the equation of motion, 
but are the intrinsic and invariant property of a Hermi- 
tian quadratic form that is equipped with commutator 
or Poisson bracket (Sec. IVIIip . 

Finally, it is worth noting that we shall confine our 
interest in this review only to the diagonalization of 
quadratic Hamiltonians. It will go out of our consid- 
eration as to whether and how a quadratic Hamiltonian 
can be derived and obtained for a real system. 



II. DIAGONALIZATION THEORY OF BOSE 
SYSTEMS 

In this section, we employ the Heisenberg equation of 
motion to study the quadratic Hamiltonian of bososns. 
It is found that a whole theory of diagonalization can be 
developed for the Bose system. 



A. Dynamic matrix 

The Heisenberg equation of motion can be derived from 
Eqs. p^ll -lLill. 



i—c 
dt 



ac + 7c^, 



. U- !■ 4. ]■ 

i—c' = — 7'c — ac' . 
at 



(2.1) 
(2.2) 



Here and hereafter, we shall apply the natural units of 
measurement, i.e., h = c ^ 1, for convenience. The two 
equations above can be combined as 



i—^jj = Dip, 
dt 



where 



D 



a 7 

— 7^ —a 



(2.4) 



It should be pointed out that the matrix D is distinct 
from the matrix M of Eq. ([TT 



M = 



a 7 
7^ a 



(2.5) 



except both a = and 7 = 0, the trivial case that has 
been excluded, ab initio, in Sec. II Al As indicated by Eq. 
p.6p . the matrix M represents the constant coefficients 
of the Hamiltonian. In contrast, Eq. (|2.3p demonstrates 
that the matrix D will control the dynamic behavior of 
the system. 

Definition 2 We shall call the matrix present in the 
Hamiltonian, such as M, the coefficient matrix, and the 
matrix present in the Heisenberg equation, such as D, the 
dynamic matrix. 



It is a characteristic feature of the Bose system that 
the dynamic matrix D is different from the coefficient 
matrix M . 

On the other hand, it can be readily seen from Eqs. 
(|2.4[) and (|2.5p that the dynamic matrix D and the coef- 
ficient matrix M have the relation, 



D = I^M. 



(2.6) 



This relation will play a fundamental role in the diago- 
nalization of the Bose system. Besides, it is worth noting 
that D is generally not Hermitian whereas M is Hermi- 
tian forever. 

Following Sec. II Bl let us study the eigenvalue problem 
ofEq. (1231) . 



ujij} = Dip. 



(2.7) 



We expect that it would generate a linear transformation 
that could be used to diagonalize the quadratic Hamilto- 
nian of bosons. 

Now that D ^ M , there arises a question. As shown 
in Eqs. (|1.20p and (|1.22p . the Hamiltonian requires a 
Hermitian congruence transformation to diagonalize the 
coefficient matrix M . However, the Heisenberg equation 
can, at most, generate a similarity transformation to di- 
agonalize the dynamic matrix D, as can be seen from Eq. 
(|2.7p . Not only the matrices to be diagonalized but also 
the manners of diagonalization are different from each 
other. That is the key problem occurring in the Bose 
system. 

To solve the problem, let us begin with a survey on the 
general properties of the dynamic matrix D. 

Lemma 3 If uj is an eigenvalue of the dynamic matrix 
D, then ~lu* will also be an eigenvalue of D. 



(2.3) Proof. The characteristic equation of Eq. (|2.7p is 



det(w/+ - D) 



Lul — a 



Lul + a 



= 0. 



(2., 



First, let us perform some elementary row and column 
operations to the characteristic determinant. 



det(a;/+ - D) 



-I 
-/ 



/ 
/ 



ujI — a —7 

7^ ojI + a 

—ujI — a — 7T 

7 —Lol + a 



/ 
/ 



(2.9) 



And then take complex conjugate. 



det(tj/+ - D) 



—uj*I — a 



—uj*I + a 



(2.10) 



where we have used the facts for the Bose system: = a 
and 7 = 7. Paying attention to 



det(w/+ - D)* = det(cj/+ - Z?) = 0, 



(2.11) 
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we have 



-u}*I — a 



= 0. 



(2.12) 



One reaches the lemma immediately by comparing this 
equation with Eq. (|2.8p . ■ 

This lemma shows that the eigenvalues of the dynamic 
matrix D will appear in pairs if they exist. When one of 
a pair is w, the other is —u*. 

Physically, this property of the dynamic matrix D orig- 
inates from the Hermitian symmetry of the Hamiltonian: 
= H . This symmetry implies that, if 

c{t) = coexp(=Fiwt) (2-13) 
is a solution of Eq. (|2.ip . then 

= 4 exp(±icj*t) (2.14) 
will be the solution of Eq. (|2.2p . That is to say, if 

7/;(f) = ■0oexp(-iwt) (2.15) 
is a solution of Eq. (|2.3p . the 

V-li) = i/'oexp(iw*t) (2.16) 
must also be a solution of Eq. (12. 3p . 



Lemma 4 // f (w) is an eigenvector belonging to the 
eigenvalue uj of the dynamic matrix D, then v{—u}*) will 
he an eigenvector belonging to the eigenvalue —uj* . Here 
the vector v{—uj*) is defined by 



with 



v{-uj*) = Y.xV*{u!) 



/ 
/ 



(2.17) 



(2.18) 



Proof. Substituting v{uj) into Eq. (|2.7p . one has 

{ujI+ - D)v{uj) = 0. (2.19) 

It follows that 

-S],(u;/+-D)S,E,w(w) = 0, (2.20) 

and that 

[-S,(c^/+ - i?)S,]* [S,i;(c^)]* - 0. (2.21) 
From Eqs. (|2.9p and (|2.10p . one can easily see that 

[-Y..,{loI+-D)Y..,]* = -LO*I+-D, (2.22) 
he thus gets 

{-LO*I+ - D)Y.^v* (w) = 0. (2.23) 



This means that T,xV*{lS) is an eigenvector belonging to 
the eigenvalue —uj* . That is just Eq. (|2.17p . ■ 

This lemma shows that, for a given pair of eigenval- 
ues {uj, —uj*), their eigenvectors can be formed into pairs 
according to Eq. (|2.17p . 

In the classical case, the dynamic matrix Dd of Eq. 
(OSD is 



D,i = KV. 



(2.24) 



It indicates that Dd is the production of a positive defi- 
nite matrix K and a Hermitian matrix V . Although Dd 
is not Hermitian in general, it is diagonalizable and all its 
eigenvalues are real. Mathematically, that is because the 
matrix K has Cholesky decomposition, as has been seen 
in Sec. II Bl Now, as shown in Eq. (|2.6p . the dynamic 
matrix D for a Bose system is the production of an indef- 
inite matrix /_ and a Hermitian matrix M . There exists 
no Cholesky decomposition for the matrix /_ , and there 
is no guarantee for D to be diagonalizable. Furthermore, 
the eigenvalues of D will, in general, be complex other 
than real even if D is diagonalizable. In a word, the 
present situation is much more involved than the classi- 
cal case. To be clear, let us take a look at the simplest 
case, i.e., the Hamiltonian of Eq. (|l.ip with n = 1. 
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H = ac) c + —7 c^c^ + —7* cc. 



(2.25) 



Solution 6 Apparently, the dynamic matrix D is a2x2 
matrix, 



D ^ 

The eigenvalue equation is 



a 7 
-7* —a 



a 7 
'7* —a 



(2.26) 



(2.27) 



Obviously, the eigenvalues can be obtained from the char- 
acteristic equation. 



uj'^-a^ + |7|^ = 0, 



the results 



±\/«'-|7l', H>l7l 
Lu^{Q, |a| = |7| 

\a\ < |7|. 



(2.28) 



(2.29) 



±n/l7l^ - a2. 



Namely, there are two real eigenvalues when \a\ > I7I, a 
zero eigenvalue when \a\ = I7I, and two imaginary eigen- 
values when \a\ < |7|. 

It is easy to show that, if \a\ — there exists only 
one eigenvector. 



(2.30) 



X 




1 


y 
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where 9 = arg(7) is the argument of 7, and the signs 
=p correspond to a = respectively. Therefore, the 

dynamic matrix D can not be diagonalized when \a\ = \ 

When I a I < I7I, the dynamic matrix D has two linearly 
independent eigenvectors. It is thus diagonalizable, but 
its eigenvalues are both imaginary. 

When \a\ > I7I, the dynamic matrix D has two lin- 
early independent eigenvectors, it is also diagonalizable. 
In particular, its eigenvalues are both real. 

In sum, the dynamic matrix of this Bose system has the 
same property as that of the classical system only when 
\a\ > I7I; It is diagonalizable, and its eigenvalues are 
real. 

This simple example exhibits clearly the complexity of 
Bose systems. To resolve this complexity, we shall study 
first the necessary and then the sufficient condition for 
the diagonalization of a quadratic Hamiltonian of bosons, 
which constitute the themes of the following two subsec- 
tions, respectively. 



B. Necessary condition for diagonalization 

To be clear and definite in the following, we would first 
give three definitions here. 

Definition 7 // an invertible matrix has the form as Eq. 
il.l9\) . we call it a Bogoliubov- Valatin matrix. 

Definition 8 A linear transformation defined by Eq. 
hi. 18]) will be called a Bogoliubov- Valatin transformation 
if T is a BV matrix and -0 is a standard bosonic field, 
i.e., it satisfies Eq. il.lS\) . 

Definition 9 The quadratic Hamiltonian defined in Eq. 
il.l\) is said to be Bogoliubov-Valatinianly diagonalizable 
if there is such a BV transformation that can fulfill both 
the conditions of Eqs. \1.21]) and \1.22]) . 

It should be pointed out that a BV transformation is 
not required to satisfy either the condition of Eq. (|f .2ip 
or the condition of Eq. (|I.22|) , according to the definition 

Lemma 10 The inverse of a BV matrix is also a BV 
matrix. 



Proof. From Eq. (|f .19p , it is easy to show that a matrix 
T is a BV matrix if and only if 



^^jfT ^-^x — T. 
Taking the inverses of the two sides, we obtain 



[T- 



This demonstrates that T -'^ is also a BV matrix. 
From Eq. (|1.7p . one can easily see that 



(2.31) 

(2.32) 
■ 

(2.33) 



This is a basic symmetry of the field operator, we shall 
call it the involution symmetry, for convenience. Mathe- 
matically, this symmetry roots from the fact that Ci and 
are not independent, but are the Hermitian conjugates 
of each other. Conversely, if a field has the involution 
symmetry as above, its component operators can not be 
independent, there must exist some relationship among 
them. 

Now, consider the new field defined by the BV trans- 
formation of Eq. (jl.lSp . It can be given by the inverse 
transformation. 



Obviously, 
which results in 

This implies that 



(2.34) 



(2.35) 



(2.36) 



(2.37) 



where Eqs. (|2.32|l - (l2.34p have been used. Therefore, the 
new field will have the same involution symmetry as the 
old one after a BV transformation. 

Lemma 11 The involution symmetry of the field opera- 
tor is conserved for the Bose system after a BV transfor- 
mation. Namely, the involution symmetry is an invariant 
property of the B V transformation. 

Suppose that ip is a new field, it thus has the involution 
symmetry. As mentioned above, its component operators 
will not be independent. In fact, it is easy to show that tp 
must have the same form as the old field defined in Eq. 
p.7p . In other words, it can be represented as follows. 



d 



(2.38) 



where 



d = 



d2 



d^^[di 4, ••• , 4 



(2.39) 



Here and c?| represent a new pair of operators, which 
are Hermitianly conjugate to each other. Their commu- 
tation rules can be derived from the inverse transforma- 
tion, 



(2.40) 
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As stated in the definition [H the old field tp satisfies the 
standard commutation rule of Eq. (|1.12|) . We therefore 
obtain 



-nt 



(2.41) 



The commutation rule for the new field ip is determined 
wholly by the BV matrix T, it may not be standard, 



(2.42) 



It is standard if and only if T satisfies the condition of 
Eq. (|1.2ip . which is equivalent to 



(2.43) 



In a word, if a BV transformation satisfies the condition 
of Eq. p.2ip . the new field is a standard bosonic field; if 
it further satisfies the condition of Eq. (|1.22p , the Hamil- 
tonian of Eq. (jl.ip gets BV diagonalized. 

We shall leave it to the next subsection to discuss how 
to obtain a BV transformation and make it fulfill the two 
conditions of Eqs. (|1.2ip and p.22p . Here and now, we 
would, above all, show a basic property of the BV trans- 
formation, which will play the central role in the theory 
of diagonalization of quantum quadratic Hamiltonians. 

Lemma 12 Under a BV transformation, the two dy- 
namic matrices respectively for the old and new fields will 
be similar to each other. 

Proof. Obviously, under a BV transformation as given 
in Eqs. (jl.lSp and (|1.19p . the equation of motion of the 
new field ip will still be linear in ip itself, i.e., 



(2.44) 



where Di is the dynamic matrix for the new field ip. Ap- 
parently, this equation has the same form as that for the 
old field Ip, 



at 



(2.45) 



where D is the dynamic matrix for the old field ip. On 
the other hand, it follows from Eq. (|1.18p that 



d , d 



(2.46) 



With the above two equations of motion for ip and (p, this 
equation can be expressed as 

D^j ^ TDiif. (2.47) 

Substituting -0 further with Eq. (jl.lSp . one has 

{T-^DT- Di)ip = 0. (2.48) 

This equation is equivalent to 

Vi(p = 0, Vie {1,2,- •• ,2n}, (2.49) 



where Vi are the row vectors of the matrix T ^DT — Di. 
It implies that 



That is 



V, =0, Vi e {1,2,--- ,2n}. 



T-^DT- Di = 0, 



Di = T-^DT. 



(2.50) 
(2.51) 
(2.52) 



In other words, the dynamic matrix will vary in a similar 
manner under a BV transformation. ■ 

Proposition 13 // a quadratic Hamiltonian of bosons 
can be BV diagonalized, then its dynamic matrix is diag- 
onalizable, and all the eigenvalues of the dynamic matrix 
will be real. 

Proof. Suppose that the diagonalized form of the Hamil- 
tonian is 



(2.53) 



i=l 



where C is a real constant. Since H is Hermitian, all LDi 
must be real, i.e., LDi G M. By definition, the new field ip 
satisfies the standard commutation rule. 



ip ■ ip' 



(2.54) 



From the two equations above, the dynamic matrix Di 
for the new field ip can be found as 



Di = diag([Ji, W2, • • • , Wn, -uJi, -W2, 



-iOn), (2.55) 



where diag(ai,a2, • • • ,a„i) denotes the diagonal matrix 
with ai, 02, • • • , Om on the main diagonal. From this, 
one can reach the proposition by the lemma [T^ ■ 

Definition 14 A dynamic matrix is said to be physically 
diagonalizable if it is diagonalizable, and all its eigenval- 
ues are real. 

By this definition, the proposition can be restated as 
follows. 

If a quadratic Hamiltonian of bosons can be BV diago- 
nalized, its dynamic matrix is physically diagonalizable. 

Corollary 15 The quadratic Hamiltonian of bosons de- 
fined in Eq. can not be BV diagonalized if the co- 
efficient submatrix a vanishes identically, i.e., a = 0. 



Proof. In this case, the dynamic matrix of Eq. (|2.4p 

reduces to 



D = 



7 
-7^ 



(2.56) 



Obviously, it is anti-Hermitian, and thus unitarily diag- 
onalizable. Since 7 7^ 0, the eigenvalues of D can not 
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all be zero, some of them must be purely imaginary. In 
other words, the dynamic matrix D is diagonalizable but 
not physically diagonalizable, which is inconsistent with 
the necessary condition for the BV diagonalization of a 
Hamiltonian. ■ 

This corollary shows that it is not all the quadratic 
Hamiltonians of bosons that can be BV diagonalized. 
What kind of Hamiltonians is BV diagonalizable? To 
answer it, one needs to study the sufficient condition for 
the BV diagonalization. 

C. Sufficient condition for diagonalization 

In the preceding subsection, we have already ob- 
tained the necessary condition for the BV diagonaliza- 
tion. Henceforth, we would presume that the necessary 
condition holds for the Bose system. Starting from this 
presumption, we shall search the sufficient condition for 
the BV diagonalization in this subsection. 

By definition, the dynamic matrix D is of size 2n. If 
D is physically diagonalizable, it has a complete set of 
totally 2n linearly independent eigenvectors. We have 
learned from the lemmas [3] and [3] that the eigenvalues and 
eigenvectors of D will appear in pairs. Let us continue 
this discussion about pairing. 

Lemma 16 // the dynamic matrix D is physically diago- 
nalizable, then, for each pair of nonzero eigenvalues, i.e., 
(lu, —uj) with Lu ^ 0, they have the same degeneracy. In 
other words, their eigenspaces have the same dimension. 

Proof. According to the lemma SI if the eigenvalue uj 
has m eigenvectors, 

vi{lu), Z = 1,2,-- - ,m, (2.57) 

the corresponding m vectors vi{—uj), 

vi{^u;)^j:,v:{io), ; = l,2,---,m, (2.58) 

are the eigenvectors belonging to the eigenvalue —uj. It 
can be readily confirmed that the vectors vi{—lu) {I — 
1,2,- • • , m) are linearly independent if and only if the 
vectors w;(w) (/ = 1, 2, • • • , to) are linearly independent. 

In other words, if to has m linearly independent eigen- 
vectors, then —to also has to linearly independent eigen- 
vectors, and vice versa. That is to say, the eigenvalues 
UJ and —UJ have the same degeneracy, their eigenspaces 
have the same dimension. ■ 

The proof above shows that, if the basis vectors for the 
eigenspace of uj {uj ^ 0) have been determined, the basis 
vectors for the eigenspace of —uj can be chosen as Eq. 
(I2.58p . and vice versa. 

Lemma 17 If the dynamic matrix D is physically diago- 
nalizable and has zero eigenvalue, the eigenspace of zero 
eigenvalue is even dimensional. In particular, its basis 
vectors can be chosen and grouped as 

v^+i{0)^J:^v^{0), Z = 1,2,---,to, (2.59) 



where 2m (m aN) is the dimension of the eigenspace of 
zero eigenvalue. 

Proof. The first point is a direct result of the lemma [TBI 
As to the second one, let us consider the eigenvalue 
equation, 

Dv{0) = 0. (2.60) 

That is a homogeneous system of 2n linear equations, its 
solution set forms the eigenspace of zero eigenvalue. 

Since the dimension of the eigenspace of zero eigen- 
value is 2to, we have rank(Z?) — 2n — 2m where rank(yl) 
denotes the rank of the matrix A. It means that the vec- 
tor v{0) has 2m free unknown components. Therefore, 
we can choose the following 2to components of v{0), 

u"(0), a = 1,2,--- ,TO,n-|-l,n + 2,--- ,n + TO, (2.61) 

as the free unknowns. First, let the free unknowns be 
respectively as follows, 

vrW^S^h (2.62) 

where 

I = 1,2, - • • ,TO, (2.63) 

a = 1, 2, • • • , TO, 71 + 1, 71 + 2, • • • , n + TO. (2.64) 

We obtain from Eq. (|2.60p the first group of eigenvectors, 

vi{0), ; = 1,2,--- ,TO. (2.65) 

Clearly, they are linearly independent. Then, using the 
lemma HI we have the other group of eigenvectors, 

v„-,+i{0) = ^,vt{0), Z = 1,2,---,TO. (2.66) 

They are also linearly independent. The definitions for 
Vm+i{0) show that 

V?n+liO) ^ S^,n+l, (2.67) 

where 

/ = 1,2, - • • ,TO, (2.68) 
a — 1, 2, • • • , TO, n -|- 1, n -|- 2, • • • , n -t- TO. (2.69) 

Equations (|2.62p and (|2.67p imply that the two groups 
are also linearly independent. The combination of the 
two groups has 2to linearly independent eigenvectors, 
they form a basis for the eigenspace of zero eigenvalue. 
All in all, the basis vectors for the eigenspace of zero 
eigenvalue can be chosen and grouped as Eq. (|2.59p . ■ 

Following the two lemmas above, if the dynamic matrix 
D is physically diagonalizable, it is enough for us to find 
a half of the eigenvectors of D, the other half can be 
determined by Eqs. (|2.58|) and (|2.59p . In other words, 
the eigenvalues and eigenvectors of D can be formed into 
pairs according to Eqs. (|2.58p and (|2.59p . Each pair 
has two linearly independent eigenvectors with opposite 
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eigenvalues. Such a pair will be called a dynamic mode 
pair. Consequently, there are totally n dynamic mode 
pairs. Henceforth, Eqs. (|2.58p and (|2.59|) will be used as 
the conventions for the pairs of dynamic modes. 

Suppose that the dynamic matrix D is physically di- 
agonalizable. One can construct a linear transformation 

I 



Here each eigenvalue is counted up to its multiplicity, and 
the n dynamic mode pairs are separated and arranged 
sequentially into the left and right halves of the matrix T. 
Since the dynamic matrix D is supposed to be physically 
diagonalizable, the matrix T has full rank and is hence 
nonsingular and invertible. This analysis demonstrates 
that an invertible linear transformation can be derived 
from the equation of motion if the dynamic matrix of the 
system is physically diagonalizable. 

I 



Paying attention to the fact that Ti^ is an elementary 
matrix, the right multiplication by it represents switching 
the left and right halves of the square matrix standing left 
to it. So we obtain 

I],T*S, = T. (2.73) 

That is an important property of the the derivative trans- 
formation, it implies that the derivative matrix T has the 
same form as that of (|1.19p . According to the definitions 
[7] and [5] as well as the lemma [TTJ we obtain the lemma 
for the derivative transformation. 

Lemma 19 // the dynamic matrix D is physically diag- 
onalizable, its derivative matrix is a BV matrix. The 
corresponding derivative transformation is a BV trans- 
formation, and thus it will conserve the involution sym- 
metry of the field operator of the Bose system. 

As shown by the proof above, it is the two conven- 
tions of Eqs. (|2.58p and (|2.59p that guarantee that the 
derivative transformation is a BV transformation. Con- 
sequently, one must comply with both of them when he 
constructs a BV transformation. 

Up to now, we have proved that a BV transformation 
can be generated by the Heisenberg equation of motion 
if the dynamic matrix of the system is physically diago- 
nalizable. 



as in Sec. II Bl 

V' = T^, (2.70) 

where represents the new field operator, and T is the 
matrix which consists of all the eigenvectors of D, 



(2.71) 



Definition 18 If the dynamic matrix D is physically di- 
agonalizable, then a linear transformation can be defined 
by Eqs. {2. 70{ l and {2. 71\j . We shall call it the derivative 
transformation, and call the corresponding matrix T the 
derivative matrix. 



By use of both the conventions of Eqs. (|2.58|1 and 
([235)1 . we have 



UJn), v(uJi), V{UJ2), v{uJn)]- (2.72) 
I 

Although the new field can inherit the involution sym- 
metry through the derivative BV transformation, its 
commutation rule is not always standard. From now 
on, we shall turn to handling this problem. As already 
known, it is determined by Eq. (|2.4ip . Usually, it does 
not matter what the magnitude of an eigenvector is. 
Above, when constructing the derivative BV transfor- 
mation of Eqs. (|2.70p and (|2.7ip . we did not consider 
the magnitudes of the eigenvectors either. Nevertheless, 
Eq. (|2.4ip shows that the magnitudes of the eigenvectors 
can change the commutation rule of the new field heavily. 
Therefore, it is necessary for us to take into account the 
magnitudes of the eigenvectors if we want to make the 
new commutation rule standard. Mathematically, the 
magnitude of a vector concerns the metric on the linear 
space. Therefore, we introduce, first, a sesquilinear form 
(l^ ITU for the Bose system. 

Definition 20 Using /_ , we define a sesquilinear form 
(p on the 2n- dimensional unitary space C^" as follows, 

\y))^{x\I^\y), y\x), \y)eC'\ (2.74) 

where Dirac notations have been used for the vectors of 
C^" . The form (j) will be directly referred to as the metric 
too. 

The definition is proper because the matrix /_ is Her- 
mitian with respect to the standard basis and standard 



T=[v{uJi), V{UJ2), v{uJn), v{~UJi), v{~UJ2), v{-UJn) ] . 

I 
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inner product of the unitary space C^". Besides, the met- 
ric vector space defined by the form (j) is nondegenerate 
because the matrix /_ is nonsingular. 

Ahhough this form is indefinite, it is very useful for 
clarifying the properties of the eigenvectors of the dy- 
namic matrix. 

Lemma 21 If the dynamic matrix D is physically diago- 
nalizable, its eigenspaces will be orthogonal to each other 
with respect to the metric /_ . 

Proof. Since the dynamic matrix D is physically diag- 
onalizable, the linear space C^" can be decomposed into 
the direct sum of the eigenspaces of D, 



= Ei®E2' 



Em, 



(2.75) 



where Ek {1 < k < m with 1 < m < 2n) are the 
eigenspaces of D, which belong to the eigenvalues tjfe G M 
respectively, 



(2.76) 



Here a label fi is added to distinguish the vectors of the 
eigenspace Ek- 

Using the relation of Eq. ()2.6p , the equation above can 
be reformulated as 

M|fc^) = . (2.77) 

As a result, we obtain 

{li^\M\kn) =uJkMI-\kn) . (2.78) 
By complex conjugate, we have 

{k^i\M\li^) =UJk{k^l\I-\l:^) , (2.79) 



the eigenvalue ujk being real. Again, from Eq. (j2.78p . we 
have 

{k^i\M\lv) =LOi{k^i\I^\liy) . (2.80) 

The combination of the two equations above leads to 

{k^i\I^\l:^) =0, ifujk^uji. (2.81) 

This equation demonstrates that the different 
eigenspaces of D are orthogonal to each other, i.e., 
Ek -L El (k ^ I), with respect to the metric /_ defined 
above. ■ 



The lemma shows that the whole space of 



■■2n 



be decomposed into an orthogonal direct sum of the 
eigenspaces of the dynamic matrix I? if Z? is physically 
diagonalizable. 



= EiQE2Q---QE„ 



(2.82) 



Meanwhile, as indicated by Eqs. (|278l) and (|2:8T|) . the 
coefficient matrix M becomes block diagonalized with 
respect to the eigenspaces of D, 



Lemma 22 // the dynamic matrix D is physically diag- 
onalizable, then, for each eigenspace of D, there exists an 
orthonormal basis with respect to the metric /_ . That is, 

{kn\ /_ \kiy) ^ X^Sf,,, \kfi) e Ek, M e Ek, (2.84) 

where 1 < ^Jl, v < nk with Uk being the dimension of the 
eigenspace Ek, and = +1 or —1. 



Proof. Taking notice of Eq. (|2.82[) . the metric /„ 
must also be a nonsingular sesquilinear form on each 
eigenspace. Otherwise, it is singular on the whole space 
C^", which leads to an evident contradiction. ■ 
Obviously, if 



Bk = {\kii) (^Ek\l<^i<nk} 



(2.85) 



is an orthonormal basis for the eigenspace Ek, then the 
union oi Bi, B2, ■ ■ ■ , B„i, i.e., 

B = 61 U 62 U • • • U S„ 

= {|fc^) (E Ek\ 1 < k < m, 1 < /i < Uk}, (2.86) 

will form an orthonormal basis for C^". 

Now, for a nonzero eigenvalue w (w 7^ 0), we can choose 
an orthonormal basis for it. 



vi{uj), 1 < I < m, 



(2.87) 



{lv\ M \kn) = ujk (kul I_ \kn) Ski- 



(2.83) 



where m is the dimension of the eigenspace of uj. Then, 
according to the lemma and the convention of Eq. 
(|2.58p . we have a basis for the eigenspace of —uj, 

vi{-uj) = J:^^v'!;{uj), l<l<m. (2.88) 

It is also an orthonormal basis, 

r "1 * 

vj{-Uj)I^Vk(-Uj) = vj{uj)Y,a;I-T,a:Vk{uj) 

= -XiSik, (2.89) 

where we have used the identity, 

E^/_S^ = -/_. (2.90) 

In sum, there are totally m dynamic mode pairs for the 
eigenenergy pair {lu, —uj) where uj ^ 0. Equation (|2.89p 
shows that each mode pair has two linearly independent 
eigenvectors with opposite norms. 

As pointed out after the lemma[Tni the two conventions 
of Eqs. (|2.58p and (|2.59p must be obeyed in constructing 
a derivative BV transformation. The discussions above 
demonstrate that for a pair of nonzero eigenvalues, the 
orthonormaliztion is compatible with the convention of 
Eq. (|2.58p . If zero is an eigenvalue of the dynamic ma- 
trix, is the orthonormaliztion compatible with the con- 
vention of Eq. (|2.59p ? i.e., does there exist such a basis 
for the eigenspace of zero eigenvalue that is orthonormal 
and satisfies the convention of Eq. (I2.59P simultaneously? 
The answer is yes. 
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Lemma 23 // the dynamic matrix D is physically diag- 
onalizable and has zero eigenvalue, there exists such an 
orthonormal basis for the eigenspace of zero eigenvalue 
that meets the requirement of Eq. \2.59\) . 

Proof. According to the lemma [171 there always exists 
such a basis for the eigenspace Vo of zero eigenvalue that 
satisfies the requirement of Eq. (|2.59p . 



Wm+fe(0) = S2;Wfc(0), /c = 1,2, • • • ,m, 



(2.91) 



where 2m (m G N) is the dimension of Vq, i.e., dim(Vo) — 
2m. 

When m = 1, 

vl{0)I.vi{0) + 0, (2.92) 
i.e., i'i(O) can not be isotropic. Otherwise, one has 

u|(0)/_wi(0) = 0, ^;t(o)/_^;2(o) ^ 0. (2.93) 
In addition, 

v\W-vM = «I(0)/-E,<(0) 
= «l(0)Ji;i*(0) 



0, 



where J is the unit symplectic matrix, 
J = 



/ 

-/ 



(2.94) 



(2.95) 



That is to say, the two eigenvectors i'i(O) and U2(0) are 
orthogonal to each other, 



?;I(0)/_W2(0) = 0, t;t(0)/_wi(0) = 0. 



(2.96) 



Equations p. 9311 and (I2.96|) show that the matrix of the 
sesquilinear form /_ vanishes identically on Vq. That 
is in contradiction with the fact that J_ is a nonsingular 
metric on Vq. In a word, the norm of v\ (0) can not vanish. 
So we can normalize it, 

u|(0)/_wi(0) = 1 or - 1. (2.97) 

Accordingly, 

wJ(0)/_W2(0) = -1 or 1. (2.98) 
Together with Eq. ()2.96|) . one has 

v|(0)/_vj(0) = -A,(5y, A, =±1, z,j = 1,2. (2.99) 

This implies that the lemma holds when to = 1. 

Suppose that the lemma holds when m = / (/ e N) , we 
shall show that it also holds when m = Z + 1. 

If there exists at least one of the eigenvectors that is 
nonisotropic, assume without loss of generality that wi(0) 
is such a vector, i.e.. 



then we have 



«t(0)/_t;i(0) ^0, 
i;/+2(0)/_«,+2(0)^0. 



(2.100) 
(2.101) 



namely, w;+2(0) is nonisotropic, too. That is because 

vi+M - S,«i*(0). (2.102) 

Consider the two-dimensional subspace W spanned by 
the linearly independent set {t;i(0), t;;+2(0)}, i.e.. 



W^ = span(wi(0),f,+2(0)). 



(2.103) 



Analogous to the case of m = 1, we can obtain such an 
orthonormal basis for W , 



4(0)/_w,(0) = -A,5y, 



(2.104) 



where Ai = ±1, and i, j = 1, / + 2. It satisfies the require- 
ment of Eq. ((2T59)) . 

If all the eigenvectors are isotropic, there must exist 
at least two eigenvectors such that their inner product is 
nonvanishing. Otherwise, the form /_ will vanish identi- 
cally on Vo, 

4(0)/_«j(0) = 0, Vi,j e {1,2,--- ,2/-|-2}, (2.105) 

which is obviously impossible because the form /_ is non- 
singular on Vo- Without loss of generality, let us suppose 
that 



z;I(0)/_«2(0) ^ 0. 



(2.106) 



We can always adjust the phase of wi(0) or W2(0) so that 
'y|(0)/_W2(0) is purely imaginary. 



z;I(0)/_z;2(0) i 



Now, set 



W2 



Vi{0)+iV2{0), 
Vl{0)-iV2{0), 



and 



Wl+2 = T,xW* = Vl+2i0) ~ ivi+3{0), 
Wi+3 = SajWj = V;+2(0) + W;+3(0). 



(2.107) 



(2.108) 
(2.109) 



(2.110) 
(2.111) 



From Eqs. (|2.106p and (|2.107|) . we obtain 

wll-wi = i\vl{0)I_V2{0)-vl{0)I-Vi{0) 

^ 0. (2.112) 

Obviously, 

wi{0), W2{0), wi+2{0), and u;,+3(0) £ Vq. (2.113) 

Besides, they are linearly independent. For convenience, 
let us reset 

vi{0) = wi, V2{0) ^W2, (2.114) 

Vl+2{0) = Wi+2, Wi+3(0) = Wi+3, (2.115) 

and consider the new set. 
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{^;i(0),i;2(0),- 



I 



,^^2i+2(0)}. 



(2.116) 



It is evident that this set forms a new basis for Vq, and 
satisfies the requirement of Eq. (|2.59p . In particular, 
vi(0) is nonisotropic, 



vl{Q)I.vi{Q) ^ 0. 



(2.117) 



Therefore, the new basis returns to the case discussed just 
above. All in all, we can always obtain a two-dimensional 
subspace W as given in Eqs. (|2.102|) - (|2.104p whether the 
basis vectors of Vq are isotropic or not. 

Using the two basis vectors of W, we can put 



C»(0) = v,+i{0) 



vl{0)I^vi{0) 



MO) 



6+,(0) = Vi+,+2{0)~ 



vl^i0)l-vi+2i0) 



«;+2(0), (2.118) 

«i(0) 



vli0)l-vi{0) 



vl^iO)I-vi+2iO) 



vi+2{0), (2.119) 



where i — 1,2,-- - ,1. They are all orthogonal to wi(0) 
and Wi+2(0), 



^;|(0)/_ej(0) = 0, i = l,Z + 2; j 



1,2,- 



,21. (2.120) 



Evidently, all those vectors ^i(O) are still linearly inde- 
pendent and eigenvectors of zero eigenvalue, i.e.. 



,21. 



(2.121) 



UO)eVo, i = l,2,. 

Now, consider the space W , 

W = span({6;(0)| i = l,2,-- - , 21}). (2.122) 

It is a proper subspace of Vq, dim(VF') = 21. Obviously, 

Vq^W®W'. (2.123) 

As shown above, wi(0) and u;+2(0) are both orthogonal 
to the set {^i(0)| i = 1, 2, • • • , 21}, therefore. 



(2.124) 



This implies that Vq is the orthogonal direct sum of W 
and W, 



It is easy to show that 



(2.125) 



i;I(0)/_«,+i(0)] ^-v\^^{Q)I^vi+,+2{0), (2.126) 
2(0)^-«.+i(0)l * = -v\{Q)I.vi+,+2{0)- (2.127) 



As a result, we obtain 

6+^(0) -Exer(O), i^l,2,--- ,1. 



(2.128) 



This indicates that the basis for W' is exacly in ac- 
cordance with the convention of Eq. (|2.59p . Since 
dim(T/F') = 21, by the induction hypothesis, the space 
W' has an orthonormal basis that satisfies the conven- 
tion of Eq. (|2.59p . Suppose the basis is the set: 



which satisfies 



where j = 1, 2, ■ 



{Ci(0),C2(0),-- - ,C2i(0)}, (2.129) 

0+,;(0)-S,Cr(0), (2.130) 
■ ■ ,1, and 

C|(0)/_0(0) = -A,<5.„ (2.131) 



where \i — ±1 and i,j — 1,2,- 
the set 



,21. It is evident that 



{z;i(0),^;i+2(0),Ci(0),C2(0),--- ,C2/(0)} (2.132) 
is a basis for Vo- Upon ordering them as follows. 



one has 



where i 



77i(0) - t;i(0), 

77;+2(0) - t;,+2(0), 

r?,+i(0) = 0(0), * = 1,2,--. 

?7/+.+2(0) = 0+.(0), i = l,'^,- 



r;,+,(0) = I],77*(0), 
• • ,1 + 1, and 
r7l(0)/_ry,(0) = -AA 



.1. 



1,2, 



(2.133) 
(2.134) 
(2.135) 
(2.136) 



(2.137) 



(2.138) 



where = ±1 and i, j = 1, 2, • • • , 2{l + 1). Here Eqs. 
^JU^, dSini, (12123, (ElSni), and ((TT?T|) have been 
used. The two equations above show that the new set 



{7yi(0),7/2(0),--- ,?72(/+i)(0)} 



(2.139) 



is an orthonormal basis for Vq , and it satisfies the require- 
ment of Eq. (|2.59p . In other words, the lemma holds for 
m^l + 1. 

Finally, by mathematical induction, the lemma is valid 
for any m E N. ■ 

In this proof, we present a modified version of 
the Gram-Schmidt orthogonalization process, which can 
maintain the satisfiability of the convention of Eq. (|2.59p . 
Put it another way, if one starts form a basis satisfying 
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Eq. (|2.59p , he can arrive finally at an orthonormal basis 
which will still satisfy Eq. (|2.59p through this modified 
Gram-Schmidt orthogonalization process. 

This lemma shows that there are totally m mode pairs 
for the zero eigenvalue w = 0, each mode pair has two 
linearly independent eigenvectors with opposite norms. 

The total n dynamic mode pairs are thus selected, their 
eigenvectors form an orthonormal basis for the whole 
space of C^", 

v\uJi)I.v{LU,)^\A3: (2.140) 

where \i — ±1 and 1 <i,j < 2n. Here each eigenvalue 
is counted up to its multiplicity. Each dynamic mode 
pair has two linearly independent eigenvectors with both 
opposite eigenenergies and opposite norms. Therefore, a 
half of the basis vectors have the norm 1, the other half 
have the norm —1. 



where, for each mode pair of (v{LL!i),v{~uJi)), the eigen- 
vectors are ordered as follows, 

v^iu;^)I^v{u;^) = 1, (2.144) 
v\~uj^)I-v{-Lu,) = -1. (2.145) 

That is to say, the left half of r„ is filled with the eigen- 
vectors with the positive norms of +1; the right half of T„ 
is filled with the eigenvectors with the negative norms of 
— 1. For convenience, we shall call Tn the normal deriva- 
tive BV matrix and Eq. (|2.142p the normal derivative 
BV transformation. To sum up, a normal derivative BV 
transformation can always be generated by the Heisen- 
berg equation of motion if the dynamic matrix of the 
system is physically diagonalizable. 

According to the stipulation for the normal BV matrix 
Tn, one has 

v''{uj,)I^v{ujj) = XAj, l<i,j < 2n, (2.146) 

where 

(1, 1 < i < n 
^^-(-1, n+l<^<2n. ^^.147) 

In terms of matrix, it can be expressed as 

T^J-Tn = I-. (2.148) 

Thereby, we arrive at the following lemma. 

Lemma 25 // the dynamic matrix D is physically diag- 
onalizable, the normal derivative BV matrix Tn satisfies 
the identity, 

Tll.Tn^I-. (2.149) 



From Eq. (I2.83P and the equation above, we obtain 
the following lemma. 

Lemma 24 With the orthonormal basis of C^" chosen 
as above, the coefficient matrix M is diagonalized, 

v\uJi)Mv{ujj) ^ \iUJi6ij, l<i,j<2n, (2.141) 

where ^^(ciJi) is an eigenvector of the eigenvalue uji, with 
Ai = 1 or —1 being the corresponding norm. 

For each mode pair, its two eigenvectors have opposite 
norms, one is -fl, the other is —1. So we can stipulate 
an order for every mode pair: The first eigenvector has 
the norm of -1-1, and the second one has the norm of — 1. 
Under this stipulation, the derivative BV transformation 
of Eqs. (P?7(I)1 and (P?7T|) becomes 



(2.142) 
(2.143) 

I 

In other words, T„ is a member of the group U {n, n) 
[l5| . The sesquilinear form /_ remains invariant under 
the transformation of Tn ■ 

This lemma implies that the new field 93 is a standard 
bosonic field, i.e., 

ip-ip^^I.. (2.150) 

Put it another way, the new component operators, di 
and dl {i — 1,2, ■■■ ,n) will satisfy the standard commu- 
tation rules for the annihilation and creation operators 
of bosons, 

[d„d]]=d,^j, [d„dj]^0, [4,4] = (2.151) 

Lemma 26 // the dynamic matrix D is physically di- 
agonalizable, the normal derivative BV matrix Tn will 
diagonalize the coefficient matrix M in the manner of 
Hermitian congruence. That is 

T^AfT„ = diag(cJi,-- - ,c^„,c^i,--- ,ujn). (2.152) 

Proof. It comes simply from the lemma and Eq. 
dSIlI]). ■ 

This lemma shows that the dynamic matrix D and the 
coefficient matrix M can be diagonalized simultaneously 
if D is physically diagonalizable: 

T^^DT„ = diag(u;i, • • • ,cj„,-wi,--- ,-uJn), (2.153) 

T^MT„ = diag(o^i,-- - ,c^„,c^i,--- ,ujn). (2.154) 

It is evident that the manners of diagonalization are dif- 
ferent: The former is diagonalized by a similar transfor- 
mation, and the latter by a Hermitian congruence trans- 
formation. In brief, the two different matrices have been 



I 

V' = Tn(p, 

Tn^[v{uJi), V(UJ2), v{uJn), v{~UJi), v{~UJ2), v{-UJn) ] , 
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diagonalized by the two different manners of transforma- 
tion simultaneously, that solves the key problem occur- 
ring in the Bose system. 

With the two lemmas above, we obtain the sufficient 
condition for the diagonalization of the Bose system. 

Proposition 27 A quadratic Hamiltonian of bosons is 
B V diagonalizahle if its dynamic matrix is physically di- 
agonalizable. 



Proof. By replacing tp with the ip of Eq. (|2.142p . Eq. 
(11.61) can be reformulated as 



(2.155) 



Using the lemma I26[ we have 

if = ^Lj,d,|d, + -^c^,--tr(a), (2.156) 



i=l 



i=l 



where, as shown in Eq. p.lSip . the di and d\ (i = 
1,2,- •• , n) are the new annihilation and creation op- 
erators for the bosons. Equation (|2.156p shows that a 
quadratic Hamiltonian of bosons can be BV diagonal- 
ized if its dynamic matrix is physically diagonalizahle. 
■ 

So far, we have found that the Heisenberg equation of 
motion is a natural generator of the normal BV trans- 
formation. This transformation can react on the Hamil- 
tonian itself and brings it into a diagonalized form auto- 
matically. That is just what we expected. 

Corollary 28 The quadratic Hamiltonian of bosons of 
Eq. (QTip is B V diagonalizahle if the coefficient suhmatrix 
7 vanishes identically, i.e., 7 — 0. 



Proof. When 7 
reduces to 



0, the dynamic matrix of Eq. 



D 



a 
-5 



(2.157) 



Since a is Hermitian, the a is Hermitian, 5^ = a, too. 
As a result, D — D'^ . It implies that D is unitarily di- 
agonalizahle and all its eigenvalues are real when 7 = 0. 
Of course, the dynamic matrix D is BV diagonalizahle if 
7 = 0, which proves the corollary. ■ 

On one hand, the corollary [^5] shows that some 
quadratic Hamiltonians of bosons are BV diagonalizahle. 
On the other hand, the corollary [15] indicates that there 
are also some quadratic Hamiltonians of bosons that can 
not be BV diagonalized. What is simultaneously the 
necessary and sufficient condition for the BV diagonal- 
ization? Evidently, the answer is just the combination of 
the two prepositions [T51 and [?fl 

Theorem 29 A quadratic Hamiltonian of bosons is BV 
diagonalizahle if and only if its dynamic matrix is physi- 
cally diagonalizahle. 



This theorem is obviously consistent with our physi- 
cal intuition: A system will behavior as a collective of 
quasi-particles if and only if there exists a complete set 
of linearly independent normal modes of motion in the 
system. 

In particular, it converts the BV diagonalization into 
the eigenvalue problem of the dynamic matrix, whose 
theory is very clear and simple mathematically, and fa- 
miliar to all of us. Therefore, this theorem makes it easy 
for us to find BV transformation and realize BV diago- 
nalization. 

Now, let us return to the example \5\ If \a\ < 
the dynamic matrix has two imaginary eigenvalues; the 
Hamiltonian is not BV diagonalizahle. If \a\ — the 
dynamic matrix is itself not diagonalizahle; the Hamilto- 
nian is not BV diagonalizahle. If \a\ > the dynamic 
matrix has two real eigenvalues, it is hence physically di- 
agonalizahle, therefore, the Hamiltonian is BV diagonal- 
izahle. To sum up, the Hamiltonian is BV diagonalizahle 
only when \a\ > |7|. 

Further, let us find out the normal derivative BV trans- 
formation for the example [5] in the case of |q;| > |7|. 

Set 



(2.158) 



The normalized eigenvector for lo can be obtained from 
Eq. (12271), 



w(w) 



1 



with the norm being 



7 

Lu — a 
7 

u! — a 



a > 
a < 0, 



(2.159) 



.*(.)/_.(.) = I ,<0. (2.160) 

According to Eq. ()2.143p . the normal BV matrix is 

T„ = (M'^)'y(-i' (2.161) 

where the normalized eigenvector for —lu can be given 
according to the convention of Eq. (|2.58p . In detail, it 
reads, 



^y2uJ{ot — UJ) 



7 to — a 
LO — a 7* 
LO — a 7 

7* LO — a 



a> 



a < 0. 



The following results can be readily verified, 
T T — J 



T^^DT — 



diag(a;, — w), a > 
diag(— w, Lo), a < 0, 



(2.162) 



(2.163) 



(2.164) 
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t1mt„ 



diag(w, ll>), a > 
diag(— w, -co), a < 0, 



_ J ^^d)d + ^uj — ^a, a > 
'Ujd^d — iuj — ha, a < 0, 



where 



[d,d^] = 1. 



(2.165) 



(2.166) 



(2.167) 



In this example, the eigenvalues are not degenerate. 
Let us look at a degenerate case. 

Example 30 

H = a{c\ci + 4c2) + jc[cl + 7*ciC2. (2.168) 
Solution 31 The dynamic matrix D is a 4 x 4 matrix, 



D 



a 7 

a 7 

-7* ~a 

-7* -a 



(2.169) 



The characteristic equation is 



2 2 I 



= 0. 



(2.170) 



The solutions are 



±^a^-\l\\ |«|>|7| 
0, \a\ = |7| 



(2.171) 



±n/l7l - a , |a| < l7l 



// \a\ < I7I, the dynamic matrix has imaginary eigen- 
values. Of course, the Hamiltonian is not BV diagonal- 
izable. 

It can be readily verified that D has only two linearly 
independent eigenvectors if \a\ — \ 



(2.172) 



where — arg(7), and =p correspond to a = ± \"f\ re- 
spectively. It means that D is itself not diagonalizable. 
Needless to say, the Hamiltonian is not BV diagonaliz- 
able when \a\ = \'y\ . 

V Iq^I > I7I; there are a pair of real eigenvalues, i.e., 
{lu,—uj) where 





- 1 " 




■ ■ 


vi{0) = 






, 1^2(0) = 


Te''» 
1 




i7r 



(2.173) 



It is found that uj has two linearly independent eigenvec- 
tors, they can be chosen and orthonormalized as follows. 



vi{uj) = < 



V2{uj) — < 






UJ — OL 

■\j2uj{a—uj) 
7 

■\J 2uj{uj — ol) 




UJ — Oi 

-y^2u;{cij — q) 





■\J 2cj(a;+a) 
■\J 2cj(a;+Q) 





UJ-\-OL 
■\J — 1uj{ijJ'\-Oi) 

7" 

yj — 2a;(ciJ+a) 





a > 



a < 0, 



(2.174) 



a > 



(2.175) 



, a < 0, 



their norms are 

vl{uj)I^vi{uj) 

{lo)I_V2{lo) 



1, a > 
-1, a < 0, 



1, a > 
-1, a < 0. 



(2.176) 
(2.177) 



The orthonormal eigenvectors for —lo, vi{—lj) and 
V2{—^), can be obtained according to the convention of 

Eq. ^Km) . 

Evidently, each eigenvalue is two-fold degenerate. The 
four eigenvectors, fi(w), V2{'jj), vi{—uj) and V2{^uj), 
form an orthonormal basis for 'C'^ , i.e., 

Tll^Tn = I-, (2.178) 

where T„ is the normal derivative BV matrix, 

[vi{uj),V2[uj),vi[-uj),V2{-uj)\, a>0 



Those facts demonstrate that the dynamic matrix D is 
physically diagonalizable if \a\ > |7|. 

To sum up, the Hamiltonian of Eq. \2.168\) is BV 
diagonalizable only when \a\ > \ 

When \a\ > I7I, the following results can be verified 
straightforwardly. 



rp-jjj, ^ r diag(u;,w, -w, ~w), a > 
" " 1^ diag(-w, -w, w,u;), a < 0, 

diag(cj, cj, w, w), a > 



diag(— — w, — w, — cj), a < 0, 



, (2.180) 
(2.181) 
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H 

where 



Ljj{d\di + c?2d2) -\- LjJ — a, a > 
—uj{d\di + d\d2) — uj — a, a < 0. 



(2.182) 



[d,,d]]^S^, 



[d„d,]=0, [dld]]^0. (2.183) 



The two examples above do not have zero eigenvalue. 
The following is an example with zero eigenvalue. 



Example 32 

H = c\ci + c\c2 — c\c2 — c\ci. 
Solution 33 The dynamic matrix D is 



(2.184) 



D = 



1-10 
-110 
0-11 
1-1 



It has three eigenvalues, 

c^ = 0, 2, 



(2.185) 



(2.186) 



For the pair of the nonzero eigenvalues, (2,-2), they 
are nondegenerate, 



1 

-1 








1 

-1 



(2.187) 



where the convention of Eq. \2. 58\) has been used. Their 
norms are 



v\2)I^v{2) ^ 1, i;t(-2)/_w(-2) 



-1. 



(2.188) 



For the zero eigenvalue, it is two-fold degenerate, the 
orthonormal basis for this eigenspace can be chosen as 
follows. 



v{0) 



1 



, ^;(-0) 



1 



(2.189) 



where the convention of Eq. h2. 59i) has been used. Their 
norms are 

v\Q)I-v{0) = 1, ^■f(-0)/_w(-0) = -1. (2.190) 

Therefore, the dynamic matrix D has three real eigen- 
values and four linearly independent eigenvectors, so it 
is physically diagonalizable. According to the theorem \29l 
the Hamiltonian of Eq. \2.184-^ is BV diagonalizable. 

According to Eq. \2.14S^ , the normal derivative BV 
matrix T„ has the form. 



T„. 



[v{2), v{0), v{-2), v{-0) 



(2.191) 



It is easy to show that 

T^I-T,,^I^, 

T-OT„ =diag(2, 0,-2,-0), 

r^Afr„ =diag(2, 0,2,0), 

H = 2d\di+0dld2, 



where 



[dud]] 



[d,,d,]=0, [did]] 



0. 



(2.192) 
(2.193) 
(2.194) 
(2.195) 

(2.196) 



For this example, the eigenspace of the zero eigenvalue 
is two dimensional, it is the smallest and nondegenerate. 
Finally, we give an example whose eigenspace of the zero 
eigenvalue is more than two dimensional. 



Example 34 



H ~ 2c|ci + C2C2 + cjc3 + V2(c|c2 + cjci) 



+ V2{c\c3 + C3C1) + (cjc3 + C3C2). 

Solution 35 The dynamic matrix D is 



(2.197) 



D = 



2 V2 V2 
V2 1 1 
V2 1 1 










-2 

-V2 

'V2 



There are three eigenvalues, 
OJ = 0, 4, 







-V2 -V2 
-1 -1 
-1 -1 



4. 



(2.198) 



(2.199) 



For the pair of the nonzero eigenvalues, (4, —4), they are 
nondegenerate. 



«(4) 



v2 
? 

! 
2 








t;(-4) 







? 



(2.200) 



with the norms as follows 

?;t(4)/_i;(4) = 1, 'yt(-4)/_w(-4) = -1. 



(2.201) 



For the zero eigenvalue, it is four-fold degenerate, the 
orthonormal basis for the eigenspace can be chosen as 
follows. 



/l 




r 1 -1 


V 3 






[2 




2\/3 


V 3 




, «2(0)= 


%/3 
2 
















. . 








(2.202) 
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vi{-0) 








1 

± 



(2.203) 



their 



norms are 



vl{0)l-vi{0) = 1, (2.204) 

wJ(0)/_W2(0) = 1, (2.205) 

u|(-0)/_wi(-0) = -1, (2.206) 

u|(-0)/_W2(-0) = -1. (2.207) 

In a word, the dynamic matrix D has three real eigen- 
values and six linearly independent eigenvectors. The 
Hamiltonian is BV diagonalizable. 

The normal derivative BV matrix T„ can he con- 
structed according to Eq. \2.14-S^ , 

T„ = [«(4),i;i(0),z;2(0),«(-4),i;i(-0),z;2(-0)] . (2.208) 

It can he verified straightforwardly that 

T^I-Tn^I-, 



T-DT,, = diag(4, 0, 0, -4, -0, -0), 
T,lA/T„ =diag(4,0,0,4,0,0), 
H = AdXdi + 04^2 + 04^3, 



whe 



[dt,dj\ 



0, [d\,d 



(2.209) 
(2.210) 
(2.211) 
(2.212) 

(2.213) 



D. Uniqueness 



The theorem [29] asserts that a quadratic Hamihonian 
of bosons can be BV diagonahzed if its dynamic matrix 
is physically diagonalizable. Apparently, there may exist 
many different BV transformations for a certain Hamilto- 
nian that can all realize the diagonalization. This occurs 
especially when some eigenvalues of the dynamic matrix 
are degenerate because there are much orthonormal bases 
for the eigenspace of a degenerate eigenvalue, and accord- 
ingly there are much various choices for the column vec- 
tors of the normal derivative BV matrix. For instance, 
it is easy to show that, when a > 0, the following two 
vectors. 



V2{uj) 



2^2w(a)+a) 
2^2w(a)+a) 



2^2uj(a-uj) 

2^2bj(u)+a) 
7* 

2^2uj(u)+a) 

gj — g 
2^2uj(a~ui) 



(2.214) 



also constitutes an orthonormal basis for the eigenspace 
of uj of the example [301 Substituting them into the BV 
matrix, one will find that this normal derivative BV ma- 
trix can diagonalize the Hamiltonian as the T„ given in 
Eq. (|2.179[) . Even if all the eigenvalues are nondegener- 
ate, the eigenvectors can choose their phases freely. In 
sum, the normal derivative BV transformation can never 
be unique. 

That poses a natural question: Can different normal 
BV transformations give rise to different diagonahzed 
forms for a certain quadratic Hamiltonian? Or, put it an- 
other way, is the diagonahzed form of a quadratic Hamil- 
tonian unique? The answer will be yes if one does not 
care the order of the quadratic terms present in a diago- 
nal Hamiltonian. 

Theorem 36 // a quadratic Hamiltonian of bosons is 
BV diagonalizable, its diagonahzed form will be unique 
up to a permutation of the quadratic terms. 

Proof. Suppose that there are two diagonahzed forms 
for the Hamiltonian of Eq. (|l.ip . According to Eq. 
(I2.156p . they can be written as 



2^^»-^tr(a), (2.216) 

i=l 



1 " 

i=l i=l 
n n 

H2 = Y1 ^l4d^ + 2 E - 2 (2-217) 

i=l 



We shall prove that the set {uji,uj2, ■ ■ ■ , w„} is identical 
to the set {lu'i,uj'2t' ' ' Here, if uji = lu'^, they are 

both counted up to the same multiplicity. 

According to the lemma I12[ the dynamic matrices of 
Hi and H2 are both similar to that of the Hamiltonian 
H of Eq. (ini). 



r-iOT„2 = D2, 



(2.218) 
(2.219) 



where Di, D2, and D are, respectively, the dynamic ma- 
trices for Hi, H2 and H, and r„i and Tn2 are both the 
normal derivative BV matrices, 



^i2^-^"2 



(2.220) 
(2.221) 



Thereby, Di and D2 are similar to each other, 

{T-^Tr.2y' Di (r-/T„2) = D2. (2.222) 

Observe 



Di = diag(wi,' 
D2 = diag(wj,' 



-Wi, • • • 
-W2,--- 



(2.223) 
(2.224) 



(2.215) We have 
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{Wl, W2, • • ■ -Wl, -W2, • • ■ , -^n} = {t^i, ^2, 



5 '^ni "'^iJ — ^^2: 



(2.225) 



where each eigenenergy is counted up to its multipUcity. 
As a consequence, 

wi e (2.226) 

or 

ujie{~uj[,^uj'^,--- ,-uj'J. (2.227) 

If 

Luie{-u;[,~cj'^,--- ,~^'n}, (2.228) 

let LOi be — w']^, i.e., uji — —Lij'i, without loss of generality. 
Thus, one obtains from Eq. (|2.222p 

v{coi) = T-^^T„2v{-io[). (2.229) 

This gives rise to 

wt(tji)/_w(wi) = v''{-u;[)I.v{-uj[), (2.230) 

where Eqs. (j2.220p and (|2.22ip have been used. However, 

v\uJi)I^v{uJi) = 1, (2.231) 
v'<{-u;[)I^v{-uj[) = -1. (2.232) 

Obviously, they contradict the equation (|2.230p . There- 
fore, 

iUii{^co[,~io'2,--- (2.233) 

it must belong to the set {uj'i,uj2, ■ ■ ■ i.e., 

toi€{co[,Lu'2,--- (2.234) 

This implies that 

{wi,W2,--- ,i^n} C {uj[,uj'2,- ■ ■ ,^'„}- (2.235) 
For the same reason, 

{^1,^2,- •• C {uji,uj2,--- (2.236) 

So 

{wi,W2,--- ,cj„} = {tji,a;2,--- (2.237) 

This means that Hi and H2 are identical up to a per- 
mutation of the members of the set {uji,uj2, ■ ■ ■ ,Wri}. In 
other words, the diagonalized form of a quadratic Hamil- 
tonian is unique up to a permutation of the quadratic 
terms. ■ 

Up to now, the BV diagonalization is always meant 
to diagonalize a Hamiltonian with respect to the nor- 
mal bosons, which fulfill the standard commutation rules. 



Sometimes, e.g., in quantum electrodynamics, the so- 
called time-polarized bosons are needed, they satisfy the 
abnormal commutation relations, 

[5„5]] = -<5y-, [6„6,]=0, [b\,b]] = Q. (2.238) 

If one exchanges the roles of bi and b\, and interprets 
them respectively as a creator and annihilator, i.e., 

d,^b\, dl^b,, (2.239) 

he has 

[d,,d]]^S,j, [d„d,]^0, [dld]] = 0. (2.240) 

That is to say, the time-polarized bosons can always 
be transformed into the normal bosons, and vice versa. 
Hence, a Hamiltonian can also be diagonalized with re- 
spect to the time-polarized bosons. Accordingly, Eq. 
(|2.156p will become 

H = - 2 - 2*'(")' (2-241) 

i=l 1=1 

where rii = —blbi are the particle- number operators for 
the time-polarized bosons. Actually, a Hamiltonian can 
be diagonalized with respect to the normal bosons, or the 
time-polarized bosons, or the mixture of both the normal 
and time-polarized bosons as you wish, e.g., 

m n 

H = '^ujidldi+ ^ iOiblbi 

i—l i—m+1 

+ 2^^^- 2 ^ (2-242) 

where < m < n. This fact will be used in Sec. IVIIIBI 
Anyway, the diagonalization is unique with respect to 
the normal bosons. The normal bosons will be used, by 
default, for the diagonalization of the Bose system unless 
otherwise specified. 

By the way, we would like to note that the conclu- 
sions up to now are also valid for the Bose system whose 
Hamiltonian is represented quadratically with regard to 
the time-polarized bosons, or the mixture of both the 
normal and time-polarized bosons. That is because, as 
mentioned above, all the time-polarized bosons can be 
transformed into the normal bosons. 

To conclude, a quadratic Hamiltonian of bosons has 
BV diagonalization if and only if its dynamic matrix is 
physically diagonalizable. If the diagonalization exists, 
its form is unique. 

Thus far, a whole theory of diagonalization has been 
achieved for the Bose system. 
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III. DIAGONALIZATION THEORY OF FERMI 
SYSTEMS 

In this section, we turn to the Fermi case. We shall 
study first the existence and then the uniqueness of the 
BV diagonalization for the Fermi system. 



A. Existence 

The Heisenberg equation for the fermionic field ijj can 
be derived from Eq. (jl.ip . 



i^V^ = Dip, 



(3.1) 



where D is the dynamic matrix for the Fermi system, 



D = 



a 7 
t -7i 



1 



(3.2) 



In contrast to the Bose system where the dynamic ma- 
trix is distinct from the coefficient matrix, the dynamic 
matrix D is now identical to the coefficient matrix M , 



M = 



7 



7 



(3.3) 



This demonstrates that the coefficient matrix M will con- 
trol the dynamic behavior of the system, just as the dy- 
namic matrix D. That is the radical difference between 
the Fermi and Bose systems. For the latter, as we know, 
the coefficient matrix does not control the dynamic be- 
havior of the system. Now that D — M and M is Her- 
mitian, D is Hermitian, too. That is another feature of 
the Fermi system, which will bring us much convenience. 

Similar to Eq. (|2.6p . the relation between D and M 
can be formally written as 



D = I.M. 



(3.4) 



This relation is useful in the diagonalization of the Fermi 
system. 

As before, let us consider the eigenvalue problem. 



tuijj = Dip. 



(3.5) 



Lemma 37 For a quadratic Hamiltonian of fermions, its 
dynamic matrix is always B V diagonalizable. 

Proof. As mentioned above, the dynamic matrix for a 
Fermi system is Hermitian. It is well known that a Her- 
mitian matrix is diagonalizable, and all its eigenvalues 
are real. Therefore, the dynamic matrix for a Fermi sys- 
tem is always BV diagonalizable. ■ 

This property is basically different from the Bose sys- 
tem. There, the dynamic matrix is not always diagonaliz- 
able. Needless to say, it is not always BV diagonalizable. 

As D = M and both are Hermitian, they can be di- 
agonalized by an exactly identical unitary transforma- 
tion. Mathematically, a unitary transformation is always 



a similar transformation, the diagonalization manner of 
D is not inharmonious with that of M any longer. The 
problem present in the Bose system disappears sponta- 
neously in the Fermi system. 

Analogous to the Bose system, one can easily show 
that the lemmas [21 HI [IHl [H and [TH are all valid for 
the Fermi system. Since the dynamic matrix is always 
BV diagonalizable now, the necessary condition for the 
BV diagonalization will hold automatically for a Fermi 
system. That guarantees further that the lemmas [iSl [13 
and [THI also hold for the Fermi system. All those eight 
lemmas stem from the Hermiticity of the Hamiltonian, 
and are irrespective of the statistics and metric of the 
system. 

By introducing a sesquilinear form with and substi- 
tuting /_ with , the lemmas [H] and [22] are both valid 
for the Fermi system, the only difference lies in that the 
norm with respect to /+ is positive definite whereas the 
norm with respect to /_ is indefinite. The lemma also 
holds for the Fermi case, but its proof needs quite a lot 
of modification, which we give below. 
Proof. The eigenspace Vq of zero eigenvalue is even di- 
mensional, let the dimension be 2m (m S N). According 
to the lemma [iTl there always exists a basis for Vq that 
satisfies the requirement of Eq. I|2.59p . i.e., 

vra+i{0)^j:^v*{0), Z = l,2,---,m, (3.6) 

When m — 1, dim(Vb) = 2, there are two basis vectors, 
i.e., wi(0) and V2{0), they are hnearly indepentdent. 
First of all, we would make wi(0) normalized, 

vl{0)I+vm = 1. (3.7) 

And then we consider 112(0). In fact, it is also normalized, 

4(0)/+z;2(0) = l, (3.8) 

that is because 

i;2(0) = S.fi*(0). (3.9) 

There are two possible cases for W2(0): (1) It is orthogonal 
to vi{0). (2) It is not orthogonal to wi(0). 
If «2(0)±«i(0), i.e.. 



we have 



vl{0)I+V2{Q) = 0, 


(3.10) 


vl{0)I+v,iO) = 1, 


(3.11) 


vl{0)I+V2iO) = 1, 


(3.12) 


vl{0)I+V2{0) = 0. 


(3.13) 


V2i0) = S.«i*(0), 


(3.14) 



and 



So the basis {wi(0), 112(0)} is itself orthonormal and sat- 
isfies the requirement of Eq. (|2.59p . In other words, the 
lemma holds if U2(0)-Lui(0). 
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If ^2(0) is not orthogonal to wi(0), i.e., 

vliO)I+V2{0) ^ 0, (3.15) 

we shall introduce two vectors wi(0) and ^2(0) as follows, 

wi{0) = avi{0) + bv2{0), (3.16) 
W2{0) = E,<(0). (3.17) 

Here a G C and b G C are two coefficients, they will be 
determined by the orthonormal conditions, 

wliO)I+wiiO) = 1, (3.18) 
wl{0)I+W2{0) = 0. (3.19) 



With this solution, we have 



Equations (|3.16p and (|3.17p show that both wi(0) and 
102(0) are linear combinations of ui(0) and ^2(0). Con- 
sequently, wi(0) and ^2(0) are also the eigenvectors of 
zero eigenvalue, i.e., wi{0) € Vq and W2{0) € Vq. 
Observe 

vl{0)I+V2iO) = [«i(0)/+I],t.i(0)]* . (3.20) 
We can adjust the phase of wi(0) anew so that 

vl{0)I+V2{0) > 0. (3.21) 
By use of Cauchy inequality, we obtain 

v\iO)I+V2{0) < ^vl{0)I+vi{0)^Jvl{0)I+V2{0), (3.22) 

where we have used the fact that vi{0) and ^2(0) are 
linearly independent. Since 



vl{0)I+vi{0)^vl{0)I+V2{0) = 1, 

we have 

t;I(0)/+i;2(0) < 1. 
In brief, we can alway have 

0<t;l(0)/+t;2(0) < 1, 



(3.23) 



(3.24) 



(3.25) 



when V2{0) is not orthogonal to vi{Q). 

Under such choice, Eqs. (|3.18p and p.l9p become 

a*a + b*b+ {a* b + b* a) vl{0) I +V2{0) = 1, (3.26) 
{a*a* +b*b*)vl{0)I+V2iO) + 2a*b* = 0. (3.27) 

It can be readily confirmed that there exists at least the 
following real solution for the coefficients a and b, 

1 



2a/1 



2a/1 



«I(0)/+«2(0) 
1 

l^l^vl{Q)I+V2{0)' 
1 

vl{0)I+V2{0) 

1 

l^l~vl{Q)I+V2{0) 



(3.28) 



(3.29) 



w;I(0)/+u;i(0) = 1, 
wl{0)I+W2{0) = 1, 

«;I(0)/+u;2(0) = 0, 
W2{0) = E,wt(0). 



(3.30) 
(3.31) 
(3.32) 
(3.33) 



That is to say, the set {u'i(O), ^2(0)} will form an or- 
thonormal basis for Vq, and satisfy the requirement of 
Eq. (|2.59p . This implies that the lemma also holds if 
^2(0) is not orthogonal to t;i(0). 

To sum up, the lemma will always hold when m — 1. 

Suppose that the lemma holds when m = Z (/ e N). 
We consider then the case where m = I + 1. Obviously, 
it has a proper subspace W spanned by the linearly in- 
dependent set {ui(0), wi+2(0)}, i.e.. 



W = spMvi{0),vi+2m. 



(3.34) 



Taking notice of 



vi+2{0) = S,«*(0), 



(3.35) 



and following the same arguments as those for the case 
of m = 1, we can obtain an orthonormal basis for W, 



(3.36) 



where = ±1 and i, j = 1,1 + 2. It is evident that this 
basis satisfies the requirement of Eq. (|2.59p . 

The rest steps of mathematical induction are simply 
similar to those for the Bose case. The convention of 
Eq. (|2.59p can be kept by the modified Gram-Schmidt 
orthogonalization process. ■ 

The lemma [M] still holds for the Fermi system, with 
Ai = 1 for i = 1, 2, • • • , 2n. 

Since all the eigenvectors are normalized to +1 now, 
one can not use the sign of the norm to stipulate an order 
within a mode pair. Here, we shall resort to the sign of 
the eigenvalue: The first eigenvalue in a pair is positive, 
and the second one negative; it is arbitrary if both the 
eigenvalues in a pair are equal to zero. Under such stipu- 
lation, the normal derivative BV transformation has the 
form, 
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T 



r 



(3.37) 
(3.38) 



where 



LJi > 0, z = 1, 2, • • • , n. 



(3.39) 



That is to say, the left half of Tn is filled with the eigen- 
vectors with nonnegative eigenvalues; the right half of T„ 
is filled with the eigenvectors with nonpositive eigenval- 
ues. 

With Tn ordered as above, the lemma [25] holds for the 
Fermi case, 



(3.40) 



i.e., T„ is a member of the U{2n) group [ISj . This lemma 
asserts that the new filed is a standard fermionic field. 
The lemma [26] must be modified as follows, 

T^AfTn = diag(cji, • • • ,u!„, -uji, ■ ■ ■ ,-t^„). (3.41) 

That is because 

TlMTn = Tll+DX, 



~ -'-n-'- + -'-n-'-n 

— T^^DT 



(3.42) 



where T^I+Tn = 1+ has been used. 

At last, we arrive at the diagonalization theorem for 
the Fermi system. 

Theorem 38 Any quadratic Hamiltonian of fermions is 
BV diagonalizahle. 

Apparently, the diagonalized form for the Hamiltonian 



n 1^1 



(3.43) 



i=l i=l 

where all the eigenenergies are nonnegative, 

tjj > 0, i = 1,2, • • • ,n. (3.44) 

Here, it is worth emphasizing that the BV diagonal- 
ization for a Fermi system is itself of unitary diagonal- 
ization, that is because r„ is, in fact, a unitary matrix, 

All in all, the BV diagonalization for a quadratic 
Hamiltonian of fermions is much simpler than that for 
a quadratic Hamiltonian of bosons. 

Example 39 

H = a{c\ci + c\c2) + l{c\cl - C1C2). (3.45) 



Solution 40 The dynamic matrix D is 



D 



a 7 

a -7 

-7 -a 

7 -a 



(3.46) 



It has only a pair of eigenvalues, (lo, —uj) where 



= \/ oP- + 7^, 



(3.47) 



they are both two-fold degenerate. The normal B V matrix 
can be chosen as follows, 

Tn = [viiuj), V2iuj), Vi{-Uj), V2i-Uj)] 
1 



[uo ~ af H 


'1^ 






7 








UJ ~ a 





-7 


UJ — a 





UJ 


— a 


7 





UJ — a 








-7 - 



, (3.48) 



where the convention of Eq. 112. 58\) has been used for 
vi{—uj) and 1)2 (—w). The diagonalized Hamiltonian has 
the form. 



H = uj{d\di + d\d2) — UJ + a. 

Example 41 

H — ^(c{c2 -t- c\ci) + l'{c\c\ — C1C2), 

where v > Q. 

Solution 42 The dynamic matrix D is 



(3.49) 
(3.50) 



D = 



fi V- 

M -t^ 

-i/ -/i 

,1^ -fi . 



(3.51) 



There are two pairs of eigenvalues, (uji , —uji ) and 
(^2, —UJ2) where 



UJi = H + ly, UJ2= f-- 1^. 



(3.52) 



The eigenvectors v{uji) and v(uj2) can be chosen as fol- 
lows, 



v{uji) = - 



(3.53) 
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Correspondingly, w(— wi) and v{—u}2) can he obtained 



from the convention of Eq. h2. 58\) . 

The normal BV matrix and the form of the diagonal- 
ized Hamiltonian are listed as follows. 

1. Iffi > u, 



){uJl), V{UJ2), v{-UJi), V{-UJ2) 



1 



H = ujid\di + U}2d\d2 - ^ (t^i + UJ2) 

2. If ti< -V, 

Tn [w(-wi), v{~uj2), v{u;i), V{UJ2) 
H = -ujid\di - uj2d\d2 + i (wi + L02) ■ 

3- If -V <ti< V, 



H = uJid\di — u!2d\d2 — - (lui — UJ2) 



(3.54) 
(3.55) 

(3.56) 
(3.57) 

(3.58) 
(3.59) 



Theorem 43 The diagonalized form for a quadratic 
Hamiltonian of fermions is unique up to a permutation 
of the quadratic terms. 



Proof. Suppose that there are two diagonalized forms 
for the Hamihonian of Eq. (|l.ip . According to Eq. 
p.43p . they can be written as 



n 1^1 

i7i = ^Mlrf,--^w^ + -tr(a), (3.60) 

1=1 4=1 

n 1^1 



where 



i=l 



uJi > 0, i — 1,2, ■ ■ ■ ,n, 
uj[ > 0, i = 1,2, • • • ,n. 



(3.62) 
(3.63) 



B. Uniqueness 

FoUowing the same arguments as those for the Bose 
system, one can readily know that there exist much dif- 
ferent BV transformations that can all realize BV diag- 
onalization to the same Hamiltonian of fermions. Nev- 
ertheless, the diagonalized form will be unique up to a 
permutation of the quadratic terms. 



J 



Let Di and D2 be the dynamic matrices for Hi and H2, 
respectively. Obviously, they are both diagonal, 

Di = diag(cJi, • • • ,LJ„,-wi, • • • ,-cj„), (3.64) 
D2 = dia.g{uj[,- ■ ■ ,uj'„,-uj'2,- ■ ■ ,-uj'„). (3.65) 

As Di and D2 are similar to each other, one has 



{W1,W2,--- ,UJn,-UJl,-UJ2, - ■ ■ = {uj[ , Uj'2 , ■ ■ ■ , Uj'„ 



(3.66) 



r 



Paying attention to Eqs. (I3.62p and (|3.63p . he can further 
obtain 



(3.67) 



This demonstrates that the diagonalized form of a 
quadratic Hamiltonian is unique up to a permutation of 
the quadratic terms. ■ 

Since all the energies of quasiparticles are nonnega- 
tive in the Hamiltonian of Eq. (I3.43p . there will be no 



J 



quasiparticle in the ground state of the system, viz., the 
ground state is exactly the vacuum of quasiparticles. If 
one does not obey the stipulation given in Eqs. (|3.37p . 
p.38p . and p.39p . he can obtain other diagonalized forms 
for the Hamiltonian. Nevertheless, the ground states of 
the system will not be the vacuum of quasiparticles. For 
example, instead of Eqs. p.37p . p.38p . and (|3.39p . let us 
stipulate anew that 



tp = Tn(p, 

I 



(3.68) 
(3.69) 



where 



> 0, i = 1, 2, • • • , n. 



(3.70) 



It is easy to show that the Hamiltonian has a new diag 
onalized form. 
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n 1 ^ 



1 



Now, all the energies of quasiparticles are nonpositive, 
the ground state of the system will be the Fermi sea which 
is occupied fully by quasiparticles. Here, the elementary 
excitations of the system would be quasiholes rather than 
quasiparticles. Upon the particlc-hole transformation. 



ft = 4> fl =dz, i = 1,2, • • • ,n, 



(3.72) 



the new diagonalized form can be transformed into the 
old one. 



n 1 ^ 



1 



i=l 



2 + -tr(a). 

1=1 



(3.73) 



and vice versa. Hence, both are essentially equivalent. Of 
course, you can also use the mixing picture if you like, 
e-g-, 



^ m 1 ^ 1 

^2^^'+ 2 ^ + (3.74) 

i—l i—m+1 

where < m < n. The elementary excitations of the sys- 
tem include now both the quasiparticles and quasiholes. 
Usually, the hole and mixing pictures are less convenient 
than the particle picture. That is the reason why we stip- 
ulate the BV transformation as in Eqs. p.37p . (|3.38p . 
and (|3.39p . By default, the particle picture will be used 
for the diagonalization of the Fermi system unless other- 
wise specified. The diagonalization is unique in the this 
picture. 

To conclude, the BV diagonalization exists and is 
unique for every quadratic Hamiltonian of fermions. 



IV. APPLICATION TO BOSE SYSTEMS 

Now, we apply the diagonalization theory of bosons 
to real systems. We shall concentrate ourselves on the 
two typical Hamiltonians: the normal Hamiltonian and 
the pairing Hamiltonian. As a matter of fact, they are 
the prototypes of many practical models, and represent 
almost all the problems which we encounter frequently in 
practice. 



A. The normal Hamiltonian 



First, let us consider the normal Hamiltonian, 



(4.1) 



In this Hamiltonian, there are only the normal terms such 
as clcj. The abnormal terms, such as cjcj and CiCj, disap- 
pear completely. This kind of Hamiltonian has been dis- 
cussed in the corollary [28l according to it, such a Hamil- 
tonian is always BV diagonalizable. In fact, the result 
can be strengthened further as follows. 

Proposition 44 The normal Hamiltonian of Eq. J^. J[ ) 

can be BV diagonalized by the unitary transformation 
generated by the coefficient matrix a. 



Proof. The eigenvalue equation is 

Dvioj) — ujv{uj), 
where D is the dynamic matrix. 



D 



a 

-5 



Let the eigenvector v{uj) be 
v{uj) = 



x{uj) 



(4.2) 



(4.3) 



(4.4) 



where x(uj) and y(w) are the two subvectors of size n. 
The eigenvalue equation of D becomes 



a 
-a 



x{uj) 



x{bj) 



It reduces to 



ax{Lo) 
ay{u}) 



LOx{uj), 

-ujy{uj). 



(4.5) 



(4.6) 
(4.7) 



The first equation is exactly the eigenvalue equation of 
the coefficient matrix a. Since a is Hermitian, it can be 
unitarily diagonalized. 



U^aU = diag(u;i,u;2, • ■ • , w„) , 



(4.8) 



where G R (i = 1, 2, • • • , n) are the eigenvalues of a, 
and U the unitary matrix which consists of the eigenvec- 
tors of a. 



U = [x[uJi), X{UJ2), ••• , x{uJn) 



(4.9) 
(4.10) 



the x{uJi) standing for the eigenvector of the eigenvalue 
cji, respectively. Obviously, Eqs. (|4.6p and (|4.7p have the 
solutions, 



V{u!i) 



x{u}i) 




i = l,2,. 



they are all the eigenvectors of Eq. 
malized as follows. 



n, (4.11) 
3, and orthonor- 



v\uJi)I^v{u}.j) ^ 5i;j, i,j 



1,2,' 



, n, 



(4.12) 
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which comes directly from Eq. (|4.9p . According to the 
convention for v{—uji), the rest half of the eigenvectors 
of Eq. gSl) is 





X*{uJi) 



, i = 1,2, • • • ,n, 



(4.13) 



they are also orthonormalized, 

v''{-UJ^)I-V{-UJJ ) ^ -S^J, i, j = 1, 2, • • • , (4.14) 

In fact, x*{LUi) is exactly the eigenvector of a, 

ax*{uji) = -{-u!i)x*{uj.^), (4.15) 

with the eigenvalue being —uji. Substituting them into 
Eq. ()2.143p . we obtain the normal BV matrix. 



U 

u* 



It can be readily confirmed that 



T T = / 



(4.16) 



(4.17) 



T„ DTn = diag(wi, • • • , w„, -wi, • • • , -w„), (4.18) 



T>fr„ =diag(wi,--- ,t^„,wi,--- ,a;„), (4.19) 



(4.20) 



where 



[rf„rf,]=0, [4,4] =0. (4.21) 



Equation (|4.16p shows that the normal BV matrix r„ 
for the diagonalization of the Hamiltonian of (|4.ip can be 
constructed using the unitary matrix U generated by the 
coefficient matrix a. Put it another way, the quadratic 
Hamiltonian of Eq. (|4.ip can be BV diagonalized by 
the unitary transformation generated by the coefficient 
matrix a. ■ 

The proposition shows that, to diagonalize the Hamil- 
tonian of Eq. (|4.ip , one should first find the unitary ma- 
trix U from the coefficient matrix a, and then construct 
the normal BV matrix r„ according to Eq. (|4.16p . 

The above procedure is feasible but somewhat redun- 
dant, it can be further simplified. As a matter of fact, the 
BV transformation corresponding to T„ can be reduced 
to a simple unitary transformation. 



(4.22) 



which can be verified readily from the substitution of r„ 
into Eq. (|2.142p . As a consequence, we obtain 



d-d) ^I, d-d^O, dt . = 0, 



(4.23) 



where the standard relations. 



c • c' 



/, c • c = 0, ct • = 0, (4.24) 



have been used. Through the unitary transformation of 
Eq. (|4.22p . the Hamiltonian of Eq. (|4.ip can be straight- 
forwardly diagonalized as follows. 



H = c^ac = d^U^aUd = ^ ujid\dj 



(4.25) 



It is the same as Eq. (I4.20p . 

Physically, this simple version arises directly from the 
fact that the Heisenberg equation of the field c. 



i—c = ac. 
dt 



(4.26) 



does not couple with its Hermitian field . If a uni- 
tary transformation U for the field c is generated by this 
equation of motion. 



Ud, 



(4.27) 



an adjoint transformation [/^ will be yielded meanwhile 
for the Hermitian field c^, 



by the equation of motion, 
.d . 



i—c 
dt 



a. 



(4.28) 



(4.29) 



As shown by Eqs. (|48l) . ([423]) and (|425l) . they diagonal- 
ize the Hamiltonian exactly. The simple version makes it 
much easier to diagonalize the Hamiltonian of Eq. (|4.ip . 
It is unnecessary to solve the eigenvalue problem of the 
dynamic matrix D and construct the BV matrix T„, but 
sufhcient for us to find out the unitary matrix U from 
the Hermitian matrix a. In short, it reduces the eigen- 
value problem of D, which is of size 2n, to the eigenvalue 
problem of a, which is of size n. 



Example 45 

H = eic\ci + e2c\c2 -)- /i(4c2 -f c\ci). 
Solution 46 The coefficient matrix is 



M £2 



It has two eigenvalues, 



UJ2 = — 

2 



Si 



e2 + ^{ei-e2) +4/^2 



£2 - V (ei - £2) + 4/i2 



(4.30) 
(4.31) 

(4.32) 
(4.33) 
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The unitary matrix can he found as follows, 

U = [ v{uji), v{uj2) ] 
1 



2 , 2 



UJ2 — £2 

uji-ei fi 



The diagonalized Hamiltonian is 

H ~ LOid\di + L02d\d2. 

Example 47 

H = e{c\ci + c\c2 + 4^3) 

+^Jl{c\c2 + c\ci + + 4c2)- 

Solution 48 The coefficient matrix is 



a 



e 
^ e fi 
e 



It has three eigenvalues, 

uji = e, UJ2 — £ + V2fJ-, — £ — V2fi. 
The corresponding unitary matrix is 



(4.34) 



(4.35) 



(4.36) 



(4.37) 



(4.38) 



U = [ v{uji), V{UJ2), v{uj3) 
V2 1 1 

V2 -V2 

-V2 1 1 



(4.39) 



(4.40) 



can 



It is easy to show that 

H = ujid\di + uj2d\d2 + wadjds. 

By the way, we note that the examples [32] and [ 
also be diagonalized using the present method. 

B. The pairing Hamiltonian 



As shown above, the Heisenberg equation and dynamic 
matrix are reducible for a normal Hamiltonian. There 
is another reducible case, which we are going to handle 
below. 

Consider the so-called pairing Hamiltonian, 

n 

H {a^ja\aj + £^jb\bj + -^^ja\b] + l*,a.fij), (4.41) 

ij = l 

where Oi {a\) and hi {h\) are both the annihilation (cre- 
ation) operators of bosons, and 



In this Hamiltonian, the abnormal terms, such as aihj 
and aj&], appear exactly in pairs, each pair has one a- 
boson and one 6-boson, there are totally n pairs between 
a- and 6-bosons. Simply speaking, the particles of the 
system are formed perfectly into boson pairs. 

According to Eq. (|2.3p . the Heisenberg equation for 
the Hamiltonian above has the variables of a; a\^ hi, and 
h\ where i = 1,2,-- - ,n, hence it has the multiplicity of 
An. It is easy to show that the Heisenberg equation can 
be reduced to the equations of motion of the variables of 
tti and 6j (j = 1, 2, • • • , n). 



. u ,+ 



(4.43) 
(4.44) 



Clearly, that is just of multiplicity 2n. We shall utilize 
those 2n multiple equations straightforwardly to study 
the diagonalization problem of the Hamiltonian of Eq. 
(|4.4ip . It is equivalent to but will be simpler than from 
the primitive equation (|2.3p and the theorem E51 as has 
already been seen from the discussions on the normal 
Hamiltonian. In particular, it will bring us a fairly simple 
algorithm for the diagonalization of the paring Hamilto- 
nian. 

For the sake of convenience, we introduce the new op- 
erators Ci and c] as follows. 



= 6 



zl ^bi, i = 1,2, - • • ,n. 



(4.45) 



Accordingly, the commutators will be 

a-a^ =1, a • a = 0, • = 0, (4.46) 
c • = -/, c • c = 0, ct • = 0, (4.47) 
a-c = 0, a-c^ = Q, a^-c = 0, a^-c^=Q. (4.48) 

In fact, the c-particles are just the so-called time- 
polarized bosons. Using those new operators, the Hamil- 
tonian of Eq. (I4.4ip can be written as 



H = xJj^MxP - tr(e), 
where M is the coefficient matrix, 



M 



and tp the field operator. 



a 7 
7''' e 



The commutator for ijj is 



(4.49) 



(4.50) 



(4.51) 



(4.52) 



a'^ — a, = e, 7 = 7. 



(4.42) 



Here, it is worth pointing out that there is no involution 
symmetry for the field ip now, which is quite different 
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from that of Eq. (|1.7p . As can be seen latter, this prop- 
erty wiU make the diagonaHzation much easier: One need 
not ensure the involution symmetry for the new field any 
more. 

The Heisenberg equation of motion for the field ip can 
be derived from Eqs. (1449]) and (|432l) . 



where D is the dynamic matrix, 



D = 



7 



-7' 



As regards M and D, one has 

D = I^M, 



(4.53) 



(4.54) 



(4.55) 



which is identical to Eq. ()2.6|) . We note that the coeffi- 
cient matrix M is Hermitian. 

Different from the normal Hamiltonian, the pairing 
Hamiltonian is not always BV diagonalizable. 

Proposition 49 The boson pairing Hamiltonian of Eq. 
^4-41\ l is BV diagonalizable if and only if the dynamic 
matrix D is physically diagonalizable. 

Proof. The sufficiency can be proved as follows. 

First, if the dynamic matrix D is BV diagonalizable, 
then its eigenspaces will be orthogonal to each other with 
respect to the metric /_ . The proof is the same as that 
for the lemma [2T1 which can be easily seen by comparing 
Eq. with Eq. (j^ . 

Second, for every eigenspace of the dynamic matrix 
D, there exists an orthonormal basis with respect to the 
metric /_ . The proof is completely the same as that for 
the lemma [211 

Now, summing up all the orthonormal bases chosen 
as above, we obtain an orthonormal basis for the whole 
space 



^2n 



(4.56) 



where v(uji) {1 < i < 2n) are the eigenvectors with = 
±1 being the corresponding norms. It follows from Eq. 
(|i35l) that 

{uj,)Mv{LUj) = XiUiSij. (4-57) 
By introducing the matrix, 

[/= [viui), V{UJ2), ••• , V{u2n)], (4.58) 

the two equations above can be formulated as 

C/t/_C/ = diag(Ai,A2,--- ,A2„), (4.59) 
U^'MU = diag(AiWi, A2CJ2, • • • , A2„c^2«). (4.60) 

Here, it is enough for us to take only into account of 
the orthonormalization of the eigenvectors because, as 



pointed out above, the field ip has no involution symme- 
try. One need not take care of both the orthonormaliza- 
tion and involution symmetry simultaneously as before, 
particularly as in the lemma Obviously, that brings 
us much convenience. 

Applying Sylveter's law of inertia to the first equa- 
tion above, we find that, of the total 2n norms (A; with 
i = 1, 2, • • • , 2n), there must be n positive norms and n 
negative norms. Upon rearranging the order of the eigen- 
vectors, v(uji), 1 < i < 2n, within the matrix [/, the two 
equations above can be reformulated as 



J/^Mt/ = diag(wi, • • • ,w„,-w„+i, 
Now, defining a new field ip, 

we have from Eq. (|4.52p 



(4.61) 
-CJ2n). (4.62) 

(4.63) 

(4.64) 



Accordingly, the Hamiltonian of Eq. (|4.49p can be writ- 
ten as 



H = if^uHlUip-iv{e). 
If one expands ip as 



LP 



t _ 



he obtains the commuatation realtions, 
d-d'' d-d^O, d'' ■d'' = 



0, 



e ■ ~ —I, e • e = 0, • = 0, 
d • e = 0, d • el' = 0, d^ ■ e ^ 0, • 
and the corresponding Hamiltonian, 

n 

H = ^{uidldi - uJn+ieUi) ~ (s) ■ 
1=1 



0, 



(4.65) 

(4.66) 

(4.67) 
(4.68) 
(4.69) 

(4.70) 



Heeding that the e-particles are the time-polarized 
bosons, we need perform the transformation. 



fi = ej, // = e^, z = 1,2, • • • ,n. 



(4.71) 



where the /-particles return to the normal bosons. At 
last, we obtain the diagonalized Hamiltonian, 

n n 
i=l i=l 

where the operators di (d|) and fi (fj) {i — 1,2, ■■■ ,n) 
satisfy the standard commutation relations for bosons, 

d-d^=I, d-d = 0, rft.rft^o, (4.73) 

f-f^^I, f-f^O, /t-/t = 0, (4.74) 

d-f^Q, d-f^^Q, d^-f^Q, dt./t^o. (4.75) 

The proof for necessity is simply similar to that for the 
proposition 1131 ■ 
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Example 50 

H = eic\ci + £24^2 + 7(44 + C2C1), 
where 7 > 0. 

Solution 51 The dynamic matrix is 



(4.76) 



D = 



£1 7 

-7 -£2 



Obviously, its characteristic equation is 

LU^ + (£2 - £i)i^ + (7^ - £i£2) = 0. 



(4.77) 



(4.78) 



1. If \ei + £2] < 27, there are two imaginary eigenval- 
ues, 



1 

"=2 



£1 - £2 ± iV 472 - (£1 + £2)" 



(4.79) 



The dynamic matrix D is not physically diagonalizable. 
2. If |£i + £2! = 27, there is only one real eigenvalue, 



1 



(£1 - £2) ■ 



(4.80) 



It is easy to show that lj has only one eigenvector. The 
dynamic matrix D is not physically diagonalizable. 
3. If |£i + £2! > 27, there are two real eigenvalues. 



UJl 



W2 



£1 ~ £2 + \/(£i + £2) "47^ 



£1 - £2 - V (£1 + £2 



472 



, (4.81) 
. (4.82) 



The dynamic matrix D is thus physically diagonalizable. 
Meanwhile, the Hamiltonian H can be BV diagonalized, 
the corresponding transformation matrix is 



U 



[v{uJi),v{ui2)] , £l+£2>27 
[v{uj2),v{uJi)], £i+£2<-27, 



whe 



v{uJi) = 

W(W2) = 

It is easy to show 



7^ - {u}i - £1) 
1 

U^I-U = 1^, 



7 

UJl - £1 

ijJ2 + £2 

-7 



(4.83) 

(4.84) 
(4.85) 

(4.86) 



U^MU = [ diag(^i, -t^2), £1 + £2 > 27 

diag(w2, -wi), £i+£2<-27, ^ ' 



H = 



ujid\di — uj2d\d-^ 
uj2d\di — uJid\d2 



2- U2- £2, £1 + £2 > 27 

UJl - £2, £1 + £2 < -27. 

(4.i 



Obviously, this example is an extension of the example 
[5(11 Besides, one can see that the present method is much 
more convenient than that adopted by the example [30l 

Example 52 

2 2 
H^Y. + + E + «^^^)' (4.89) 

i—l 

where £ > 0, and 7 > 0. 

Solution 53 The dynamic matrix is 



D 



■ £ 7 7 " 

£ 7 7 

—7 —7 — £ 

— 7 —7 — £ 



There are four eigenvalues, 

UJl = £, tJ2 = -£, ^^3 = t^+i ^^4 = 

w/iere 

±jv^472 - £2, £ < 27 

W± = < 0, £ = 27 

±v/£2 -472, £>27. 



(4.90) 



(4.91) 



(4.92) 



If £ < 27, _D /las two imaginary eigenvalues, H is not 
BV diagonalizable. 

If £ = 27, the zero eigenvalue has only one eigenvector, 
H is not BV diagonalizable. 

If e > 2^ , D is physically diagonalizable, H can be BV 
diagonalized. The corresponding transformation matrix 
is 



U=[v{uJi), V{UJ3), V{UJ2), V{UJ4) 



(4.93) 



whe 



v{uji) = 



V{UJ2) 



1 

V2 










1 

V2 



V{UJ3) 



27 



V87^-2(a>+-£)^ 
27 

V87^-2(a)+-e)^ 
V87^-2(c^+-e)' 
V87"-2(a>+-£)^ 



(4.94) 



t;(a;4) = 





V87^-2(t^- 

iiJ- +£ 




^87^-2(1^- 
27 




V87^-2(w 
27 




\/8t'-2('^ 





The diagonalized Hamiltonian is 



(4.95) 



H 



ujidldi + ujzd\d2 — ^21^3^3 — ujid\di 
—UJ2 ~ UJ4 — 2e. 



(4.96) 



The propositions and algorithms developed in this sec- 
tion can be applied to statistical as well as condensed- 
matter physics 0, H, S i, [H, [B [3 . 
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APPLICATION TO FERMI SYSTEMS 



Solution 56 The coefficient matrix is 



As in tiie preceding subsection, we shall also concen- 
trate ourselves on the normal and pairing Hamiltonians. 
They represent the problems which we encounter most 
frequently in practice. 



A. The normal Hamiltonian 

The normal Hamiltonian reads 



e 

It has two eigenvalues, 

0^1 = £ + /i, L02 = e — ^. 
The unitary matrix can be easily found, 



U ■ 



1 



1 1 

1 -1 



4 ,j = l 



(5.1) 



1. If e > fi, LOi > and L02 > 0, 



present case. 

Proposition 54 A normal Hamiltonian of fermions can 
he BV diagonalized by the unitary transformation gener- 
ated by its coefficient matrix. 

That is also because the Heisenberg equation for the 
Hamiltonian of Eq. (|5.ip is reducible. It can be reduced 

as 





Cl 


1 


'l 1 ■ 




'di' 


easily transplanted to the 


C2 


~ 71 


1 -1 







i—c = ac. 
dt 



(5.2) 



Since a is Hermitian, this equation of motion can gener- 
ate a unitary transformation U for the field c, 



where 



c = Ud, (5.3) 

C/tj7 = UU^ = I, (5.4) 
U^aU = disLg{uJi,uj2,- ■ ■ ,uj„). (5.5) 



The d represents the new field, it is easy to show that d 
is a standard fermionic field, 

d-d^^I, d-d^O, d^-d^^O. (5.6) 



Accordingly, 



(5.7) 



Besides, a particle-hole transformation will be needed if 
some eigenenergies are negative. 

To sum up, a normal Hamiltonian can be BV diagonal- 
ized by the unitary transformation generated by its co- 
efficient matrix no matter whether the system is bosonic 
or fermionic. 

Example 55 

H = e{c\ci -I- C2C2) -I- ti{c\c2 + c\ci), (5.8) 
where > 0. 



H = LiJid\di + u}2d\d2. 
2. If —fi < e < ^, uji > and 002 < 0, 



Cl 


1 


■ 1 1 ■ 










1 -1 




.4. 


H = 


ujid\di 


- L02d\d2 


+ ijJ2 



3. If £ < -fi, < and LlI2 < 0, 



Cl 


1 


■ 1 1 ■ 




\d\] 


C2 


"71 


1 -1 




.4. 



4. ± 

H = —LOid[di — aJ2C?2^2 + CJi + LJ2- 



(5.9) 
(5.10) 
(5.11) 

(5.12) 
(5.13) 

(5.14) 
(5.15) 

(5.16) 
(5.17) 



Here, a particle-hole transformation is performed to the 
d2- particles if — /i < e < 11, and to both the di- and d2- 
particles if £ < —fJ.- 



B. The pairing Hamiltonian 

The pairing Hamiltonian reads 

n 

H ^ ("y"4ai + + lijatbj + 7j>!^^)' (5-18) 



where 



= a, = e, 7 = —7. 



(5.19) 



Following the Bose case, let us introduce the new op- 
erators Ci and cl as 

^t J 



Ci = b\, c] = bi, i = 1,2, - • • ,n. 



(5.20) 

The new anticomniutators will be 

a-a^^I, a-a = 0, 0^-0^=0, (5.21) 
c-c^^I, c-c = 0, ct.ct = 0, (5.22) 
a • c = 0, a • = 0, • c = 0, • = 0. (5.23) 
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Obviously, they are still standard, which is rather differ- 
ent from the Bose case. In terms of these new operators, 
Eq. (|5.18p can be expressed as 



H = tp'^ Mip + iv{e), 
where M is the coefficient matrix, 



M 



a 7 
7^ — £ 



and if) the field operator, 



It is evident that 



, [at, ct]. 

M = M\ 



(5.24) 



(5.25) 



(5.26) 



(5.27) 
(5.28) 



Namely, M is Hermitian, and i/) is a standard fermionic 
field. This means that Eq. ()5.24p is, in fact, a normal 
Hamiltonian. Hence, we obtain the proposition. 

Proposition 57 A pairing Hamiltonian of fermions can 
be first transformed into a normal Hamiltonian, and then 
B V diagonalized by the unitary transformation generated 
by the corresponding coefficient matrix. 



As well known, Eq. (|5.20p represents a particle-hole 
transformation in physics. So, the proposition says ac- 
tually that a pairing Hamiltonian of fermions can be 
transformed into a normal Hamiltonian by a particle- hole 
transformation. 



Example 58 

H = eic{ci + e2c\c2 + l{c\c\ + C2C1), 
where 7 > 0. 

Solution 59 The coefficient matrix is 



(5.29) 



M 



£1 7 
7 -£2 



It has two eigenvalues 
ei - £2 + 



UJ2 



£1 + £2) 



-472 



£1 - £2 



£i+£2 



-472 



(5.30) 

(5.31) 
(5.32) 



and generates a unitary transformation matrix, 
U = [ v{uji), V{UJ2) ] 

1 



\/i^i-£iy 



+ j2 



7 LU2+ £2 

Wi - £1 7 



(5.33) 



If 1^ ^ — £i£ and £1 > 0, the diagonalized Hamilto- 
nian has the form, 



H = ujid\di + Lij2d\d2 + £2, 



where 



di 
d2 



= u- 



ci 



(5.34) 



(5.35) 



2- If 7^ < — £i£2 o,nd £1 < 0, the diagonalized Hamil- 
tonian has the form, 



H = —ujidldi — ijj2d\d2 + ijJi + ijj2 + £2, (5.36) 



where 



d\ 

4 



Cl 



(5.37) 



3. If 7^ > — £i£2, the diagonalized Hamiltonian has 
the form. 



H = ujid\di — uj2d\d2 



- W2 + £2, 



where 



di 
l4 



Cl 



(5.38) 



(5.39) 



Obviously, this example is an extension of the example 
1391 One can see that the present method is much simpler 
than that used by the example [39l 

Example 60 

2 2 
H^Y. + + E ^^^b] + a,b,), (5.40) 

i—l iii — 1 

where £ > 0, and 7 > 0. 

Solution 61 The coefficient matrix is 



M 



£07 7 ■ 

£ 7 7 

7 7 — £ 

7 7 — £ 



It has four eigenvalues, 

LOl = £, UJ2 



(5.41) 



(5.42) 
(5.43) 



= ^£2 + 472, = -^£2 + 472. 

The diagonalized Hamiltonian is 

H — ujid\di — L02d\d2 + uj^dl^d^ — u^d^di 

4-t^2 + W4 + 2£. (5.44) 

The corresponding B V transformation has the form. 
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4_ 





1 





1 


27 


V2 
1 


27 


V2 





\/87^+2(w3-e)'' 


V2 

1 


^872+2(w3-e)'' 
ti^3— e 


V2 


^872+2(w3-e)^ 



V87^+2(w4+e)^ 

\/87^+2(w4+e)^ 
27 

^872+2(^4+5)^ 
27 

^872+2(^4+6)^ 



4 



4 



(5.45) 



The propositions and algorithms developed in this sec- 
tion can be applied to statistical p hysi cs [1 , H, 0, H, 0, H] , 
condensed-matter physics [l^, US l2ll |. and nuclear 
physics 



A. Theory of Dirac diagonalization 

Consider the Hamiltonian that is quadratic in coordi- 
nates and momenta, 



VI. GENERALIZED BOGOLIUBOV-VALATIN 
TRANSFORMATION 

Historically, Dirac fH and later Shrodinger [H, [13, [H, 
[26j found that a linear harmonic oscillator can be diag- 
onalized with regard to the bosonic creation and annihi- 
lation operators. 



H 



2m 

uja^a - 



2 2 
-muj q 



(6.1) 



where m and lo are the mass and (angular) frequency of 
the oscillator, respectively. The q and p are the coordi- 
nate and momentum operators. They are both Hermi- 
tian, and satisfy the canonical commutation rule, 



[p, q] = 



(6.2) 



The a and are the annihilation and creation operators. 
They are Hermitian conjugates of each other, and satisfy 
the bosonic commutation rule, 

[a,at] = l. (6.3) 

Both a and are the linear functions of q and p, 



at = 



V 2 V mu!/ 
Imu! ( 



muj/ 



(6.4) 
(6.5) 



and vice versa. For convenience, we shall call such a 
transformation Dirac transformation, and call the corre- 
sponding diagonalization Dirac diagonalization. 

In this section, we shall study the necessary and suffi- 
cient conditions for Dirac diagonalization. We find that 
Dirac diagonalization is actually a generalization of BV 
diagonalization. 



" ^ 1111 



If If If 1 t t 

-p'^ip+ -q'Kq+ -^p'7q+ 2? 7 



(6.6) 



where the coefficients are all real, nu G M, fj,ij G M, 7ij € 
R, and 



II = ll\ K — . 



(6.7) 



The coordinates and momenta satisfy the canonical com- 
mutation relations, 

[Pt,qj\ = [Pt,Pj\ = 0, [qi.qj] = 0, (6.8) 

or equivalently, 

p-q'^^-il, p-p^^O, q-q-i^O. (6.9) 

Let us put 



The Hamiltonian can written more compactly as 
where M is the coefficient matrix, 



M 



7 

7^ K 



which is Hermitian. Equation (|6.9p becomes 
6 ■ 6^ = E„, 



where 



-il 
il 



Obviously, Y.y is a Hermitian matrix. 



(6.10) 



(6.11) 



(6.12) 



(6.13) 



(6.14) 
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The Heisenberg equation of motion for the field (f) can 
be derived from Eqs. (|6.1ip and (|6.13p . 

i^6 = D^, (6.15) 
at 

where D is the dynamic matrix, 

D = Y.yM. (6.16) 

GeneraUy, D is not Hermitian. The above equation gives 
the relationship between the dynamic matrx D and the 
coefficient matrix M, it is the counterpart of Eq. (|2.6p . 
The eigenvalue equation for Eq. (|6.15p is, 

ujv{lj) = Dv{uj). (6.17) 

It can be rewritten as 

iujv{uj) = JMv{uj), (6.18) 

where J is the unit symplectic matrix given in Eq. (|2.95p . 
Paying attention to the fact that JM is a real matrix, one 
realizes immediatly that the lemma [3] also holds now. 

Lemma 62 If lo is an eigenvalue of the dynamic matrix 
D, then ~uj* will also be an eigenvalue of D. 

Nevertheless, the lemma |4] should be modified as fol- 
lows. 

Lemma 63 // v{uj) is an eigenvector belonging to the 
eigenvalue oj of the dynamic matrix D, then its complex 
conjugate v*{uj) will be an eigenvector belonging to the 
eigenvalue . 

These two lemmas show that the dynamic mode pair 
appears now in the form of {lo,v{lj)} and {— w*, i;*(w)} 
where the two eigenvectors are complex conjugates of 
each other. 

Now, consider the Dirac transformation, 

(jj^T^, (6.19) 

where represents the bosonic field given in Eq. (|1.7p . 
Evidently, the Heisenberg equation of ijj is still linear and 
homogeneous, 

i^^ = Dii^, (6.20) 

where Di is the dynamic matrix for ■0. Following the 
proof for the lemma [T^ we obtain the lemma below. 

Lemma 64 Under a Dirac transformation, the two dy- 
namic matrices respectively for the old and new fields will 
be similar to each other. 

It is evident that Di is a real diagonal matrix if the 
Hamiltonian of Eq. (|6.1ip has been diagonalized Dira- 
cianly. As a result, the proposition [T3] holds for Dirac 
diagonalization, too. 



Proposition 65 // a Hamiltonian quadratic in coordi- 
nates and momenta can be Diracianly diagonalized, its 
dynamic matrix is physically diagonalizable. 

Also, it can be readily confirmed that the lemma [H] 
holds for the present case. 

Lemma 66 // the dynamic matrix D is physically di- 
agonalizable, then, for each pair of nonzero eigenvalues, 
i.e., {uj,—uj) with lo ^ Q, they have the same degener- 
acy. Namely, their eigenspace have the same dimension. 
Especially, their bases can be chosen as 

vi(-Lo)^v*i(u;), l = l,2,--- ,m, (6.21) 

where m (m (z N) is the dimension of the eigenspace of 

LO. 

But the lemma [T71 needs a few modifications. 

Lemma 67 // the dynamic matrix D is physically diago- 
nalizable and has zero eigenvalue, the eigenspace of zero 
eigenvalue is even dimensional. In particular, its basis 
vectors can be chosen and grouped as 

vm+i{0)^vt{0), ; = l,2,---,m, (6.22) 

where 2m (m (zN) is the dimension of the eigenspace of 
zero eigenvalue. 

Proof. It just needs to prove the second point. 
When w 0, Eq. ((O^ becomes 

JMv{0) = 0. (6.23) 

As pointed out above, JM is a real matrix, its eigen- 
vectors can be chosen naturally as real vectors. In other 
words, there exists a real basis for the eigenspace of zero 
eigenvalue, 

wi{0)^w^{0), Z = 1,2,- • • ,2to. (6.24) 

Let us put 

vi{0) = wi{0)-^iw,n+i{0), (6.25) 
vm+i{0) = wiiO)-iw^+i{0), (6.26) 

where Z = 1, 2, • • • , m. They are obviously a basis that is 
in accordance with Eq. (|6.22p . ■ 

Instead of /_ , we can here introduce a sesquilinear form 
using Ey, it is also a nonsingular metric. Notice the sim- 
ilarity between Eqs. (|6.16p and (|2.6p . We can transplant 
the lemmas [2T] and [22l to the present case. 

Lemma 68 // the dynamic matrix D is physically diago- 
nalizable, its eigenspaces will be orthogonal to each other 
with respect to the metric Sy. 

Lemma 69 // the dynamic matrix D is physically diag- 
onalizable, then, for each eigenspace of D, there exists an 
orthonormal basis with respect to the metric Yiy. 
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As to the lemma [23l it needs some modifications. 

Lemma 70 // the dynamic matrix D is physically diag- 
onalizable and has zero eigenvalue, there exists such an 
orthonormal basis for the eigenspace of zero eigenvalue 
that 

vm+i{0)=vt{0), / = l,2,---,m, (6.27) 

where 2m (m G is the dimension of the eigenspace of 
zero eigenvalue. 

Proof. According to the lemma l67l the eigenspace Vq of 
zero eigenvalue has a basis, 

Vm+i{0) = vtiO), l = l,2,--- ,m, (6.28) 

where dim(Vo) = 2m (m G N). 

When m ~ 1, wi(0) must be nonisotropic, 

vl{0)T,yVi{0) ^ 0. (6.29) 

Otherwise, one has 

wJ(0)S,yWi(0) = 0, t;t(0)Sj,V2(0) = 0. (6.30) 
In addition, 

vimyV2{o) = vimyvm 

= -ivl{0)Jvl{0) 

= 0. (6.31) 

That is to say, 

wl(0)EyW2(0) = 0, 4(0)SyUi(0) = 0. (6.32) 
I 



where 

v^uj,)^yv(uj,) = 1, (6.39) 
v\-uj^)J:yv{-uj^) = -1. (6.40) 
Consequently, we obtain a new version of the lemma 

m 

Lemma 71 // the dynamic matrix D is physically diag- 
onalizable, its derivative transformation Td satisfies the 
identity, 

Tl^yTd = I-. (6.41) 
The lemma shows that 

iP-^^^I.. (6.42) 



Equations (|6.30p and (|6.32p contradict the fact that Ej, 
is a nonsingular metric on Vq. That is to say, vi(0) can 
not be isotropic. It can thus be normalized, 

w|(0)Ej,t;i(0) = 1 or - 1. (6.33) 

Correspondingly, 

vl{0)T.yV2iO) = -I or 1. (6.34) 

Those together with Eq. (|6.32p show that 

z;J(0)Ej,t;j(0) = -A,(5y, A, = ±1, i,j-l,2. (6.35) 

Thereby, the lemma holds when m = I. 

The rest steps of mathematical induction are similar 
to the lemma [23l ■ 

The combination of the lemmas (pS|) " ([7(I| shows that 
there are totally 2n dynamic mode pairs. Each mode 
pair takes the form of {a;,u(w)} and {—lu,v*{lu)}, viz., 
it contains two opposite eigenenergies, and two complex 
conjugate eigenvectors. Most importantly, the two com- 
plex conjugate eigenvectors have opposite norms, one is 
+1, the other is —1: 

V^-UJ^)'SyV{-UJ^) = - [w^ (wi)Ej,u(u;j)] * = ±1. (6.36) 

Thereby, we can stipulate an order for every mode pair 
as in Eqs. ()2.142|) - (I2.145|) : The first eigenvector has the 
norm of +1, the second one has the norm of —1. So, we 
obtain a derivative transformation Td as follows, 



(6.37) 
(6.38) 

I 

The new field -0 has another important property. 

Lemma 72 The new field ij} difined in Eq. \6.3T^ has 

the involution symmetry, 

i,= {^^^)\ (6.43) 
Proof. As is a real field, one has 

(/) = ^. (6.44) 
From Eq. (fOTjl . it follows that 

= TfiA, (6.45) 



4> = TdV', 

W(CJ„), v{~-UJi), V{-UJ2), v{~UJn)], 
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which means 



It is easy to show that 



(6.46) 
(6.47) 



Substituting it into Eq. (|6.46p . one obtains 

= Td^x^^. (6.48) 

Since Td is invertible, and S"^ — Tt^, he arrives finally 
at 



(6.49) 



It is equivalent to Eq. (|6.43p . ■ 

This lemma shows that the new field ip takes exactly 
the form of Eq. (|1.7p . Together with Eq. (|6.42|) , it means 
that the new field is a standard bosonic field. 

Lemma 73 // the dynamic matrix D is physically diago- 
nalizable, its derivative transformation Td will diagonal- 
ize the coefficient matrix M in the manner of Hermitian 
congruence, 

T^MTd = diag(a;i, • • • , Ci;„, wi, • • • , w„). (6.50) 



Proof. From Eq. (|6.16|) . it follows that 



Therefore, 



tUit,, 



M = Y.yD. 



- Tl^yDTd 



(6.51) 



(6.52) 



= Tl^yTdT^'DTd 
= I^T^^DTd. 

As D is physically diagonalizable, we have 

^d^D^d = diag(LJi, • • • ,LJ„, -cji, • • • , -ujn)- (6.53) 

The combination of the two equations above proves the 
lemma. ■ 

Theorem 74 A Hamiltonian quadratic in coordinates 
and momenta is Diracianly diagonalizable if and only if 
its dynamic matrix is physically diagonalizable. 

Proof. The necessary condition has been proved by the 
propostion l65l 

The sufficient condition can be proved as follows. 

Consider the quadratic Hamiltonian of Eq. (|6.1ip . If 
its dynamic matrix is physically diagonalizable, it can 
generates a derivative transformation as given in Eq. 
(|6.37p . Under this transformation, Eq. (|6.1ip becomes 



H 



-ij^TlMTd^ 

n 1 ^ 



(6.54) 



where Eqs. pTT)) . ((02)) . and ((6?50)) have been used. It 
is a Hamiltonian that is diagonal with respect to the 
bosonic operators d\ and (i = 1, 2, • • • , n). ■ 

This theorem shows that the derivative transformation 
defined in Eq. (|6.37p is exactly a Dirac transformation, 
it brings a Hamiltonian quadratic in coordinates and mo- 
menta into the form diagonalized with respect to bosons. 
Just as BV transformation, Dirac transformation can be 
generated by the equation of motion of the system auto- 
matically. 

Now, let us return to the linear harmonic oscillator, 
the dynamic matrix is 



D = 



-i 

1 



i 

mLO^ 



It has a pair of eigencnergies. 

If w = 0, there exists only one eigenvector, 



^;(0) 



It means that the Hamiltonian of a free particle, 

,,2 



H = 



P 
2m' 



(6.55) 



(6.56) 



(6.57) 



(6.58) 



is not Diracianly diagonalizable. If lu > 0, Z? is physi- 
cally diagonalizable, there are two linearly independent 
orthonormalized eigenvectors 



u(-cj) = 



1 



V2mw 
1 



V2r 



—imuj 
1 

imuj 
1 



(6.59) 
(6.60) 



where the convention of Eq. (|6.2ip has been used for 
v{—Lo). They generate a Dirac transformation. 



(6.61) 



V 


1 


—imuj 


imoj 




a 


q 


y/2muj 


1 


1 




at 



As a result, we have 

P = 

q = 



-ixl—-{a - a'), 



^/2'r 



= (a + at). 



(6.62) 
(6.63) 



They are just the inverse of Eqs. (|6.4p and (|6.5p . Ac- 
cording to Eq. (|6.54p . the diagonalized Hamiltonian has 
the form. 



H = uja^a H — w. 



1 



(6.64) 



which is the same as Eq. (j6.1 
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Theorem 75 // a Hamiltonian quadratic in coordinates 
and momenta is Diracianly diagonalizable, its diagonal- 
ized form will he unique up to a permutation of the 
quadratic terms. 

Proof. The proof is similar to that for the theorem 



Set 



Ci = -j={qt + ipi), 



V2 



{q, -ipi), 



(6.65) 
(6.66) 



where qi and pi (i = 1, 2, • • • , n) are the coordinates and 
momenta which satisfy the canonical commutation rules 
of Eq. (|6.9p . It is easy to show that ci and c[ are the an- 
nihilation and creation operators that satisfy the bosonic 
commutation rules of Eqs. (|1.2p - (|1.4p . The inverse is also 
true. Upon such an invertible linear substitution, the two 
Hamiltonians of Eqs. and (|6.6p can be transformed 

into each other, up to a real constant. Thereby, the BV 
diagonalization of Eq. p.ip is mathematically equiva- 
lent to the Dirac diagonalization of Eq. (|6.6p . However, 
Dirac diagonalization does not require, in physics, the 
initial field to be a standard bosonic field. Correspond- 
ingly, the transformation matrix does not need to ensure 
the invariance of the metric, and is no longer limited to 
be a member of the group U{n, n), cf., 



TII^T„ = I-, 



(6.67) 
(6.68) 



In this sense, Dirac diagonalization is more general than 
BV diagonalization. In particular, it can be generalized 
easily to complex collective coordinates and momenta, 
so it is rather useful in the quantizations of Bose fields, 
which will be discussed in the next section. 

Here and now, we would give an interesting example 
of Dirac diagonalization. 



that the magnetic field B is set along the z-axis. Thus 



1 



1 



2m 

-UJLixPy 



ypx 



2m 



Pz 



(6.71) 



where uj^ = qB /2m is the Larmor frequency. Obviously, 
Pz commutes with iJ, it will be conserved. This means 
that Pz can be replaced by a constant, which brings a 
contantant energy, p^/2m, to the Hamiltonian. Since 
a constant energy is unimportant to a Hamiltonian, we 
shall concern ourselves with the simplified version, 



H 



1 



t^L{xPy - VPx 



1 



2(2 



y 



(6.72) 



As shown in the following, this Hamiltonian can be Dira- 
cianly diagonalized and yield the so-called Landau levels. 



Solution 76 The dynamic matrix is 



D = EyAf, 



where M is the coefficient matrix, 



M = 



The dynamic matrix D has three eigenvalues 
uji = 2wl, ^2 = ~2ujl, Us = 0. 



m 








-UJL 





m 












muj\ 













muj\ 



(6.73) 



(6.74) 



(6.75) 



The first two constitute a dynamic mode, their eigenvec- 
tors are complex conjugate to each other. 



v{2uj) 



2^/mujL 



-mujL 
-imujL 
—i 
1 



(6.76) 



B. Landau quantization 

Consider a charged particle moving in a uniform mag- 
netic field. The Hamiltonian is 



H 



(6.69) 



where q, m and p denote the charge, mass and momen- 
tum of the particle, respectively. As to A, it is the vector 
potential of the magnetic field, which can be expressed 
in the form [27l. 



A = X r, 
2 



(6.70) 



where B represents the magnetic field, and r the coordi- 
nates of the particle. Assume without loss of generality 



v{-2lu) = 



2^mujL 



—mujL 
imtOL 

i 

1 



the corresponding norms are 

V^2uJL)'^yV{2uJL) 
v\-2uOL)'i^yV{-2LOL) 



-1. 



(6.77) 



(6.78) 
(6.79) 



The third is zero, it is two-fold degenerate, and has two 
linearly independent eigenvectors, e.g.. 



vM = 



mujL 


1 



f2(0) = 





-mujL 
1 




(6.80) 
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They can be linearly combined and orthonormalized into 
a dynamic mode, 



v{+0) 



1 



-iniLOL 
i 
1 



imuL 
—i 
1 



(6.81) 



with the norms being 

v\-0)Y.yv{-0) = -1. 



(6.82) 



(6.83) 
(6.84) 



The Dirac transformation is therefore obtained as fol- 
lows, 



Px 

Pv 

X 

y 



ai 

0.2 

J 



Finally, the diagonalized Hamiltonian is 
H — 2L0La\ai + uj^. 



(6.85) 



(6.86) 



This result is rig orous and exactly the same as that 
due to Landau |27| . with lUc = 2uol being the cyclotron 
frequency. 



C. Partial diagonalization 



If 7 = 0, 



(6.87) 



it reduces to the Hamiltonian of Eq. (|1.24p . As aheady 
known from Sec. II Bl the latter is identical to the Hamil- 
tonian of Eq. (|1.56[) . Thereby, we obtain the proposition. 



Proposition 77 The Hamiltonian of Eq. |g.<g?| ) is Dira- 
cianly diagonalizable when the matrices fj, and k are both 
positive definite. If fj, is positive definite but k is nonneg- 
ative definite, it is partially diagonalizable, i.e., except 
that the normal modes with zero frequencies are conserved 
quantities and hence not Diracianly diagonalizable, all 
the rest part can be Diracianly diagonalized. 

The Hamiltonian of such kind is quite important be- 
cause it is frequently encountered in various physical 
problems, particularly, in small oscillations and field 
quantization. To handle it, one can first transform it 
into Eq. (|1.56p and then perform a Dirac diagonalization 



to the latter, the total Dirac transformation being the 
product of the two successive transformations. That will 
be easier than handling Eq. (|6.87p directly. 

In the end, it is worth noting that all the conclusions of 
this section are valid for the time-polarized commutation 
relations, 



(6.88) 



which can be seen readily by exchanging the roles of the 
group of Pi (Vi e S" C {1, 2, • • • , n}) and the group of qi 
(yi € S C {1,2,-- - ,n}). This kind of abnormal com- 
mutation relations occurs in the quantization of Maxwell 
field, and will be discussed in Sec. IVIII Bl 



VII. FIELD QUANTA 

In this section, we intend to examine field quanta, in- 
cluding Klein-Gordon field, phonon field, and Dirac field. 
In references, such problems are less addressed by BV or 
Dirac transformation. We find that they are pretty good 
tools for those problems. 



A. Klein-Gordon field 

Let us begin with the neutral Klein-Gordon field 4>{x) 
[H, HSIlOl- Its Hamiltonian reads as follows, 

H = y dxi [n'i^] + [V</)(x)]' + m2</,2(x)} , (7.1) 

where m > is the mass of the field. The 7r(x) is the mo- 
mentum density conjugate to the field 0(x), they satisfy 
the canonical commutation rules. 



[7r(x),0(x')] = -z<5(x-x'), 
[^(x),7r(x')] = 0, 
[0(x),0(x')] = 0, 



(7.2) 
(7.3) 
(7.4) 



where (5(x) denotes Dirac delta function. 

As usual, we would expand (/'(x) and 7r(x) into plane 
waves, 



(x) 



1 



dp(/)(p)e 



2p-X 



(7.5) 



TT X 



(27r)^/^ 

^ / dp^(p)e-P-, (7.6) 



i2n) 



where we used the duality between 0(x) and 7r(x). Phys- 
ically, the (j){p) and 7r(p) represent the complex collective 
coordinates and momenta of the system, respectively. As 
(/)(x) and 7r(x) are both real-valued fields. 



4>Hp) = 4>{-p), 

7rt(p) 7r(-p). 



(7.7) 
(7.8) 
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In terms of 0(p) and 7r(p) , the Hamiltonian and canonical 
commutation rules can be expressed as follows, 

H = y dpi [7rt(p)7r(p) + (m2 + p^) (/.t(p)0(p)] , 

(7.9) 

k(p),0(p')] = (7.10) 
[7r(p),7r(p')] = 0, (7.11) 
[<^(p),(/.(p')] = 0. (7.12) 

First of all, let us take a look at the equations of motion 
for 0(p) and 7r(p), which can be derived from Eq. (|7.9p 
and Eqs. (frTUl) - ([71^ . 



i7rt(p), 



i^7r1'(p) = -i(m2+p2) 



They can be combined as 
d 



*^¥'(p) = D(p)(p(p), 



where (/?(p) is the field operator, 

7rt(p) 



P) 



</'(p) 

and £>(p) the the dynamic matrix 
D{V) 



-i (m^ + p2) 

1 



Obviously, the commutation rule for (y5(p) is 

'^(p)-^^(p') -^(p-pOs, 



-'y^ 



(7.13) 
(7.14) 

(7.15) 
(7.16) 

(7.17) 
(7.18) 



where Sj, is defined in Eq. (|6.14[) . 

The form of Eq. (|7.15p hints us that the Hamiltonian 
H should be formulated with the field </'(p)j 



H= I dpi^t(p)M(p)^(p) 



where M(p) is the coefficient matrix, 

1 



M(p) = 
It is evident that 



TO^ + p^ 



^(-P) - [^^P)], 

D(p) = E^Af(p) 



(7.19) 



(7.20) 



(7.21) 
(7.22) 



With those relations on hand, the following conclusions 
become self-evident. 



1. The Hamiltonian of Eq. ()7.9p is Diracianly diago- 
nalizable if and only if D{p) is physically diagonal- 
izable. 



2. If D{p) is physically diagonalizable, the dy- 
namic mode pair takes the form of {lo,v[uj)) and 
(oj, v{—uj)) where w(— w) = v*{uj). 

3. If D{p) is physically diagonalizable, there exists an 
orthonormal basis with respect to the metric "Ey. 
It generates a Dirac transformation as follows. 



Td = [ w(w), v{-oj) 



where 



v'^ {uj)Y,yv{u}) = 1, 
w^(— [j)Sj,w(— tj) = —1. 

This implies that 

and that 

^(p)-V'^(p') ='5(P-P')^- 



(7.23) 
(7.24) 

(7.25) 
(7.26) 

(7.27) 

(7.28) 



4. The new field ■0(p) has the involution symmetry, 

t 



V'(-p) = (£^)) 



(7.29) 



This together with Eq. (|7.28p means that V'Cp) is 
a standard bosonic field, and assumes the following 
form, 



V'(P) 



a(p) 



where 



[a(p),at(p')] = ^(p-p'), 

[a(p),a(p')] = 0, 
[at(p),at(p')] = 0. 



(7.30) 



(7.31) 
(7.32) 
(7.33) 



Those discussions demonstrate that the theory of Dirac 
diagonalization is also suitable for the complex collective 
coordinates and momenta. 

It is easy to show that -D(p) has a pair of real eigen- 
values. 



±£(p) = ±-\/ m? + p2. 



(7.34) 



Since D{p) is a square matrix of size 2, this means 
that D{p) is physically diagonalizable. According to 
the above statement 1, the neutral Klein-Gordon field 
is Diracianly diagonalizable. 

According to the above statements 2 and 3, the corre- 
sponding orthonormal eigenvectors can be chosen as 



v{e{p)) = 
v{~e{p)) = 



1 



\/2i(p) 
1 



-ie[p) 
1 

ie{p) 
1 



(7.35) 
(7.36) 
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with the norms being 



B. Phonon field 



T;t(£(p))S,«(e(p)) = 1, (7.37) 
z;t(-e(p))S,«(-£(p)) = -1. (7.38) 

Finally, according to the above statements 3 and 4, the 
Dirac transformation can be constructed as follows, 



7rt(p) 



wi(£(p)), ■i;2(-£(p) 



a(p) 

at(-p) 



(7.39) 



As a result, we have 



(P) 



2£(p) 



[a(p) + at(~p)] , (7.40) 



Hp) = -^^^[a(-p)-at(p)], (7.41) 



and 



dp 



e(p)a^(p)a(p) + ie(p) 



(7.42) 



Substituting Eqs. ([7^ and into Eqs. ((731) 

and (|7.6p . and complementing them with the variable of 
time, we have 



(x,t) 



1 



'^'"V 2(2^)^e(p) 
at(p)e-»[p-x-e(p)t] 



a(p)e 



i[p-x-e(p)t] 



(7.43) 



7r(x,t) 



-i I dpJ-^M_<{a(p)e'[p-x-.(p)*] 



2(27r 
-at(p)e-'[P-''-^(P)*l 



(7.44) 



They are explicitly Lorentz covariant. 

Evidently, all those results are the same as Refs. [H, 
[2^, [sO]. The charged Klein- Gordon field can be handled 
similarly. 

Besides, we note that other complete sets of orthonor- 
mal functions, e.g., spherical waves, can be used instead 
of the plane waves to expand the fields if necessary. 

For the neutral Klein- Gordon field, one can deal with 
the real collective coordinates and momenta as in Sec. 
IVII if he expands the fields 0(x) and 7r(x) with the real 
plane waves, i.e., {sin(p-x), cos(p-x)}. He can transform 
back to the complex representation at the end of the cal- 
culation. That is rather tedious. In quantum mechanics 
and quantum field theory, complex waves and fields are 
unavoidable. That is the reason why we generalize the 
diagonalization theory given in Sec. IVII to the case of 
complex collective coordinates and momenta. 



For the neutral Klein-Gordon field, Eq. (|7.9[) shows 
that it belongs to the case that the matrices /z and k in 
Eq. ()6.87|) are both positive definite. Physically, that is 
because it is massive, i.e., to > 0. There are also fields 
that belong to the other case where fj, is positive definite 
but K is nonnegative definite. A familiar example is the 
phonon field. 

For the sake of brevity, we shall consider a simple three- 
dimensional lattice. As usual, we take the harmonic ap- 
proximation under which the Hamiltonian of 
the system becomes 

+ ^EE'J'"/3(1-1')""(V(1'), (7.45) 

l.a 

where to is the mass of the atoms or ions, a and (3 denotes 
the a;-, or y-, or z-component, 1 is shortened for the lattice 
vector. The u(l) represents the phonon field, and p(l) the 
conjugate momentum field. They satisfy the canonical 
commutation rules, 



r(i),/(i')] = 0, 

K(l),u'5(l')] = 0. 



(7.46) 
(7.47) 
(7.48) 



The matrix <I>q;3(1 — 1') stands for the interaction between 
the atoms or ions. 



$a/3(l-l') 



du^{\)dup{\') 



$;3a(l'-l), (7.49) 



where $ is the elastic potential of the lattice, and the sub- 
fix denotes the equilibrium configuration of the atoms 
or ions. Since the equilibrium configuration corresponds 
to the minimum of the potential, the matrix <I'c(/3(l — 1') 
is nonnegative definite. Moreover, it fulfils the condition. 



^<i>aMl-l') = 0, 



(7.50) 



due to the translational symmetry of the lattice. 

Using the collective coordinates u"(k) and momenta 
p"(k). 



p"(k) 



1 



ikl 



kl 



(7.51) 
(7.52) 



we obtain 



H = ^Eb"W]tp"(k) 

+\ $aMk)K(k)]+u^(k), (7.53) 



k,a,/3 
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where 







(7.54) 


b"(k),/(k')] 


= 0, 


(7.55) 




= 0, 


(7.56) 



and 



<i>„;5(k)-^$„^(l)e 



-ikl 



(7.57) 



Here and hereafter in this subsection, all the vectors k 
belong to the first Brillouin zone. 

Since a simple lattice always has inversion symmetry, 
^a/3(k) gets the following properties. 



$„Mk) = $;^(k), 

<^al3{k) = $/3a(k), 
'^>ap{-k) = $a/3(k). 



(7.58) 
(7.59) 
(7.60) 



The first two properties indicate that ^apik) is a real 
symmetric matrix. Besides, ^afsik) is also a nonnegative 
matrix, just as $q,^(1). 

The Hamiltonian of Eq. (|7.53p can be handled using 
the field. 



^(k) 



p\k) 
p^k) 

p3(k) 

K(k)]' 
[-^(k)]' 

k3(k)lt 



(7.61) 



as in the proceeding subsection. However, it will be more 
convenient to follow the propositions[T]and[77l We there- 
fore perform, first, a linear transformation that will make 
H diagonalized with regard to the new collective coordi- 
nates and momenta. As in Sec. II Bl this transformation 
can be produced by the equations of motion of the coor- 
dinates itQ,(k), exactly speaking, the following eigenvalue 
equation, 



u;^u°'(k) 



1 



-5^$„^(k)u'^(k). 



(7.62) 



As <I>Q,^(k)/m is a real, symmetric, and nonnegative ma- 
trix, it has three eigenvalues, 

w^(k)>0, cr==l,2,3, (7.63) 

and a complete set of three orthonormal eigenvectors, 

'^-"^1'^-" ^'-^ - ^ ■ (7.64) 



a. 

E 



e-(k)e^,(k) = S,^,, 
e^(k)e^(k) = Sap. 



(7.65) 



They are the so-called polarization vectors. One can fur- 
ther adjust these eigenvectors such that 



e„{k) = e^(-k). 



(7.66) 



that is because ^api^k) = ^ap{k). According to Eqs. 
([TM)) and prsS)) . u"(k) andp"(k) should be expanded 
as 

u"(k) = 5]0.(k)^e^(k), (7.67) 
p"(k) = ^^,(k)V^e^(k), (7.68) 

a 

where 4)cr{k) and Haik) are the new collective coordi- 
nates and momenta of the system. In terms of these new 
collective coordinates and momenta, Eq. (|7.53p can be 
expressed as 



H 



E 



^t(k)^^(k)-f ic.2(k)4(k)0„(k) 



where 



[7r,(k), (k') 
[^,(k),7r,,(k') 
[0,(k),0,,(k') 



and 



Uk) - 



0, 
0, 



.(-k), 4(k)=7r^(~k). 



(7.69) 



(7.70) 
(7.71) 
(7.72) 



(7.73) 



The Hamiltonian of Eq. (|7.69p is the same in form as 
that of Eq. (|7.9p . it can thus be treated as the latter. 
However, the photon field is massless. 



,(k)^0, k^O, 



(7.74) 



which can be readily seen from Eqs. (|7.50p and (|7.57p . 
Therefore, the Hamiltonian H is only partially diagonal- 
izable, the components of k = can not be diagonalized. 
For k ^ 0, the Dirac transformation is the same as Eqs. 
and (EH]), 



r(k) = 



1 



2a;,. (k) 



[a.(k)H-at(-k)], (7.75) 



TTa{k) = 



r(k) 



[a.(-k)-at(k)] . (7.76) 



Under this transformation, the Hamiltonian becomes 

a 

+ EE ^<T(k)a[^aka + ^w„(k) 



. (7.77) 



That is the partially diagonalized form for the Hamilto- 
nian of the photon field. 

It is easy to show that ^^{G) {a = 1,2,3) represent 
physically the momenta of the center of mass of the sys- 
tem. The partial diagonalization is thus not difficult to 
understand because the center of mass of the system be- 
haviors as a free particle. 
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Remark 78 Usually, the components o/ k = are re- 
garded to be diagonalizable as those o/k ^ 0. That is to 
say, the Dirac transformation. 



1 



[a,(k)+at(-k)], (7.78) 



and the Hamiltonian 



M/^[«.(-k)-at(k)], (7.79) 



H = Y^ w<,(k)a[^ak^ + iw^(k) 



k.cr 



(7.80) 



are taken to be appropriate for all k fSll . \33i . \33i l. Strictly 
speaking, that is not right, particularly in the discrete case 
(Notice that LUa-Q<.)\]i=o = 0. The transformation becomes 
meaningless when k — 0.). Nevertheless, it will cause no 
problem in the thermodynamic limit. That is because the 
state of\i — has zero measure and contributes nothing 
to the integration over k. You can change the values of 
an integrand on a null set at your will, that imposes no 
influence on the integration. Since one usually needs to 
take the thermodynamic limit finally in his calculation, it 
will be convenient to think that the total Hamiltonian is 
diagonalizable in this case. The same thing occurs in the 
photon field Isdi l and in the magnons of antiferro- 

magnets fr\ . \J8]. In the language of mathematics, those 
fields can be said to be diagonalizable almost everywhere. 

All the discussions in this subsection are also valid for 
the diagonalization of the Maxwell field under Coulomb 
gauge. The diagonalization of the Maxwell field under 
Lorentz gauge will be discussed in Sec. IVIII Bl 



C. Dirac field 

Finally, let us consider the Dirac field [H, [H, [s^l . It 
is a Fermi field, different from the two cases above. The 
Hamiltonian is 



H 



dx-!/;^(x) (— ia • V + /3m) V'(x), 



where m is the mass of the field, and 



<J^ 
<7^ 



/ 

-/ 



(7.81) 



(7.82) 



with ai [i — 1,2,3) being Pauli's 2x2 matrices, and / 
the 2x2 unit matrix. The spinor fields, i'i'x.) and ■0^(x), 
satisfy the anticommuation rules, 

[V'^(x),Vi(x')] = <5(x-x')5,,., (7.83) 
[V'^(x),V.(x')] = 0, (7.84) 
[V'i(x),V;t(x')] = 0. (7.85) 



As before, we can expand the fields V'(x) and V-'U^) 
with plane waves, 



^(x) 
V'^(x) 



(27r 



3/2 



(27r 



,3/2 



dpV(p)e''P-^, 
dpV^(p)e^'P-' 



(7.86) 
(7.87) 



where 



V'(p) - 



^t(p) 



(7.. 



■ ci(p) 

C2(P) 
C3(P) 
. C4(p) 

cl(p), 4(p), 4(p), 4(p)]. (7.89) 



Substituting them into Eq. (|7.8ip . we obtain 

i/ = y dp-0t(p)(p.^_^^^)^(p) 

= y dp ^ c|,(p) (p • a + m/3)^, c,(p), (7.90) 



where 



[c^(p),4(p')] = <5(p - p')<5^., 

[C;.(p),C,(p')] = 0, 

[ct(p),4(p')] = 0. 



(7.91) 
(7.92) 
(7.93) 



Since a. 



ai and /J^ = /?, _ff is a normal Hamilto- 
nian that has been discussed in Sec. IV AI According to 
the proposition [SH it can be diagonalized by the unitary 
transformation, 



V'(p) = Tp^(p), 
where Tp is the unitary matrix. 



(7.94) 



(7.95) 



Tt (p . a + m[3) Tp = diag(e(p), e(p), -£(p), -£(p)). 



and iy9(p) the new field, 



di(p) 

d2(p) 

4(-p) 
4(-p) 



(7.96) 



(7.97) 



Here, a particle-hole transformation has been performed 
for the negative energies. 

After the transformation, H becomes 

H= / dp^e(p)4(p)d^(p), (7.98) 
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where 



K(p),4(p')] 
[4(p),4(p')] 



^(P - P')Sf,,y, 

0, 
0. 



(7.99) 
(7.100) 
(7.101) 



As usual, the vacuum energy has been removed from the 
Hamiltonian [H, [H, (s^l . 

In this section, the theory of transformation has been 
apphed to examine totally three kinds of field quanta. As 
has been seen, it operates neatly and concisely. 



VIII. MATHEMATICAL ESSENCE OF DIAGO- 
NALIZABILITY 

Up to now, we regard BV or Dirac diagonalization as a 
physical consequence of the Heisenberg equation of mo- 
tion. Diagonalization represents the normal modes of 
motion of a system. Such a view provides us a concrete 
picture and intuitive interpretation of diagonalization, 
and thus makes it easy to understand. However, there is 
yet another view, it is abstract but more fundamental. In 
this view, the diagonalization in itself is an intrinsic and 
invariant property of a Hermitian quadratic form that is 
equipped with commutator or Poisson bracket, neither 
Hamiltonian nor equation of motion will be needed any 
more. 



A. Heisenberg operator 

For simplicity, we shall take the Hamiltonian of Eq. 
()l.ip as an instance to explain this abstract view. 

First, we regard Eq. (jl.ip purely as a Hermitian form 
Q that is quadratic in creation and annihilation opera- 
tors, 

n ^ ^ 

Q= ("y4cj + 2^v44 + ^IjiCtCj). (8.1) 

And then we introduce a linear operator / to replace the 
Heisenberg equation of motion. 



f{x) = i[Q,x]. 



•2) 



where 



= '^^{ziCi + Sic|), Zi £ C, and Si e C. (8.3) 



For convenience, we shall call / the Heisenberg operator. 
Obviously, we can define such a Heisenberg operator for 
each Hermitian quadratic form. 

If one denotes the fith component of the field -(/; defined 
in Eq. (|1.7[) with i/j^, he obtains 



where are exactly the entries of the dynamic matrix 
D. In matrix notation, it can be written as 



(8.5) 



which is the counterpart of Eq. (|2.3p . 

By replacing H with Q and the Heisenberg equation 
with the Heisenberg operator /, one can readily show 
that all the lemmas, propositions, and theorems in Sec. 
im will still hold. Evidently, it is also true for the other 
cases, e.g., the Fermi system. Here and now, there is 
no Hamiltonian, no time, and no equation of motion. 
This demonstrates unambiguously that the BV or Dirac 
diagonalization is essentially the algebraic property of a 
quadratic Hermitian form. From this standpoint, any 
physical picture and interpretation are redundant and 
unnecessary, they have nothing to do with the self of 
the BV or Dirac diagonalization and can be completely 
removed away. 

It follows immediately from this abstract view that all 
the conclusions of Sec. IVII hold for classical systems, up 
to a real constant. 

Such sublation enlarges the scope of the objects of di- 
agonalization, it can be performed to all the Hermitian 
quadratic forms besides the physical Hamiltonians. To 
show this point of view more straightforwardly, let us 
look at the following example. 
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Jz = xpy - ypx 



(8.6) 



As well known, Jz is the orbital angular moment along 
the z-direction. It is quadratic in coordinates, i.e., x 
and 2/, and momenta, i.e., px and py. Since Jz is not 
a Hamiltonian, the Heisenberg equation of motion be- 
comes meaningless. Nevertheless, it can be Diracianly 
diagonalized. 



Solution 80 The dynamic matrix is 



D = Y.yM, 



■7) 



where M is the coefficient matrix, 



M = 



-1 
10 
10 
-10 



There exists a pair of eigenvalues for D, 
OJ = ±1. 



(8.9) 



It is easy to show that each eigenvalue has two linearly 
independent eigenvectors, for example, 



(8.4) 



^-2(1) = 



(8.10) 
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They can be linearly combined and orthonormalized as 
follows, 



«i(l)=2 



(8.11) 



with the norms being 

z;t(l)E,z;i(l)-l, «|(l)I],i;2(l) = -1. 
Correspondingly, the Dirac transformation will be 



.12) 



Px 

Py 

X 

y 



= [vi{l),V2{-l),V,{^l),V2{l)] 



ai 

0.2 



(8.13) 



where 

Finally, the diagonalized form of Jz is 
Jz — a\ai — a\a2. 



(8.14) 



This result is quite similar to the coupled boson rep- 
resentation for angular momentum due to Schwinger 
[H, HI]. It is rigorous, concise, and agrees exactly with 
the familiar result about the orbital angular momentum, 
i.e., the eigenvalues of Jz consists of all the integers. 

There is another interesting example. 
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H = -- 



2m 



2 2 

-muj q , 



(8.16) 



whe 



> and LU > 0. 



Formally, it looks like a "negative harmonic oscillator" . 
Of course, it can not represent a real physical system. 
Thereby we treat it just as a Hermitian quadratic form. 
It is ready to conjecture that this quadratic form would 
be diagonalized, by the same Dirac transformation as for 
the normal harmonic oscillator. 



Solution 82 The dynamic matrix is 



D = SyM, 



where M is the coefficient matrix, 





M = 







.17) 



(8.18) 



The dynamic matrix D has a pair of eigenvalues, 

£ = ±L0. (8.19) 



and two orthonormal eigenvectors, 

1 



v{uj) 

v{-uj) 



1 



y/2r 



imuj 
1 

^imuj 
1 



v'^ {uj)Y,yv{uj) = — 1, 
v\~ijj)Y,yv{^uj) = 1. 

They generate a Dirac transformation. 



1 



= (a + a'l'). 



8.20) 
8.21) 

8.22) 
8.23) 

8.24) 
8.25) 



They are the same as Eqs. \6. 62\) and 116. 63\) . The diag- 
onalized form is 



H = —Lja^a uj. 

2 



(8.26) 



(8.15) Those results confirm our conjecture completely. 



This example reminds us that the two propositions [T] 
and [77] also hold for the case where /i is negative def- 
inite and K is nonpositive definite, no matter whether 
the commutation rules are standard, or time-polarized, 
or mixing. This fact will be useful in the following sub- 
section. 



B. Time-polarized photons 

In this subsection, we shall concern ourselves with the 
diagonalization of the Maxwell field under Lorentz gauge. 

As well known, the Maxwell field contains two spurious 
degrees of freedom in Lorentz gauge, viz., the longitudi- 
nal and time-polarized components. Because these de- 
grees of freedom are not of physical quantity, the Hamil- 
tonian is merely a Hermitian quadratic form. The ab- 
stract view of diagonalization finds its way and role here. 

In Lorentz gauge, the Hamiltonian of Maxwell field has 
the form [28, 29, 30], 

F = - i y dx [7r'^(x)^^(x) + V^^x) • V^^(x)] , 

(8.27) 

where v4^(x) {fj, — 0,1,2,3) are the vector potentials, 
and 7r^(x) the corresponding canonical conjugate fields. 
They satisfy the following commutation rules, 

[Anx),^''(x')] - *5^''<5(x-x'), (8.28) 
[^^(x),A''(x')] = 0, (8.29) 
K(x),^'^(x')] = 0, (8.30) 
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where g'^'^ is the metric tensor, 



where 



1 











-1 



.31) 



Equations (|8.28p - (|8.30p indicate that the Hamiltonian H 
belongs to the mixing case where there are both the stan- 
dard and time-polarized commutation relations simulta- 
neously. Physically, that roots from the requirement of 
Lorentz covariance. 

Expanding these fields with plane waves. 



..32) 
.33) 



we obtain 



i/ = -i I dx{K(p)]t^^(p)-Hp^ [A^ip)]U,ip)}, 

(8.34) 



where 



[^''(P),^''(P')] 
K(p),A''(p')] 
[7r^(p),^''(p')] 



and 



[A'^(p)]t 

K(p)]t 



^g'^'Sip - p'), 

0, 

0, 



^n-p): 
^n-p). 



(8.35) 
(8.36) 
(8.37) 



(8.38) 
(8.39) 



Just like the phonon field, the polarization vectors can 
be obtained fr-om the eigenvalue equation for the fields 
^'^(p), 



LO^cj^ip) = D{p)cj,{p), 



where 



^(P) 



A^ip) 
A\p) 

aHp) 



(8.40) 



,.41) 



and 



D{P) 











..42) 



Paying attention to the fact that D{p) — p^I, the or- 
thonormal basis for the polarization vectors can be cho- 
sen, except p = 0, as follows. 



e(p,0) = (1,0,0,0), 
e{p,i) = (0,ej(p)), 



(8.43) 
(8.44) 



e(p,3) = ^. 



..45) 



e,(p)-ej(p) = %, i = 1,2,3. (8.46) 

That is to say, e(p, 0) is the time-like polarization vector, 
e(p,i) (i = 1,2) the space-like transverse polarization 
vectors, and e(p, 3) the space-like longitudinal polariza- 
tion vector. Obviously, 



e^(p,A)e^(p,A') =5AA' 



..47) 



i.e., the polarization vectors form a four-dimensional or- 
thonormal system. Using this new basis, the Hamiltonian 
can be expressed as 



^ = I fdp\j2{i^>^ip)]^^>^ip)+p'iMp)]^Mp) 



ko(p)]^^o(p)+p' [Ao(p)]tAo(p)^ 



where 



and 



[^a(p),^A'(p')] 

[A^{p),Ay{p')] 

[7rA(p),7rA'(p')] 



[Mp)]" 

['rA(p)]t 



i.9AA"5(p 

0, 

0, 



^a(-p), 
7rA(-p)- 



(8.49) 
(8.50) 
(8.51) 



(8.52) 
(8.53) 



It is worth noting that the time- like components, Aq{p) 
and 7ro(p), consititute a "negative harmonic oscillator". 

According to Sec. IVII A[ H can be diagonalized by the 
following Dirac transformation. 



^a(p) 
ttaCp) 



1 



where 



and 



2£(p) 



£(P) 



aA(p) + 4(-p) 
aA(-p) -oUp) 



iPh 



[ax{p),a\,{p')] = -g\y5{p-p'), 
[ax{p),ay{p')] 0, 
[ai(p),4,(p')] = 0. 



(8.54) 
(8.55) 

(8.56) 



(8.57) 
(8.58) 
(8.59) 



Here, as usual, the commutation rules for the transverse 
and longitudinal polarizations (A — 1, 2, 3) are chosen as 
normal or standard; those for the time-like polarizations 
(A = 0) are chosen as abnormal or time-polarized. The 
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abnormal bosons here are usually called scalar or time- 
polarized photons. That is also the reason why we call 
the particles satisfying Eq. (|2.238p the time-polarized 
bosons. The diagonalized form of the Hamiltonian is 



H= / dp£(p) 



^ 4(p)«a(p) - aS(p)ao(p) 



,A=1 



(8.60) 



where the vacuum energy has been removed from the 
Hamiltonian. This equation holds in the sense that the 
Hamiltonian is Diracianly diagonalizable almost every- 
where in momentum space (except p = 0). 
Finally, the quantized fields can be written as 



A^'{x) 



dp 



1 



^ 2(2^)^£(p) 
aA(p)e'^-" + 4(p)e-'P-" 



aA(p)e^^-^-4(p)e-'f^ 



.61) 



.62) 



It can be readily seen that all those results are the same 
asRefs. [Ullill^. 

The diagonalization of the Maxwell field under Lorentz 
gauge is rather complicated. First, it contains unphys- 
ical degrees of freedom. Second, it needs mixing com- 
mutation relations, both for initial and final fields. And 
third, it is not diagonalizable everywhere but almost ev- 
erywhere. One sees that those problems can be resolved 
naturally by the diagonalization theory developed in this 
review. 



is equipped with commutator or Poisson bracket. 

The theory is dynamic matrix oriented. 

The dynamic matrix can be derived from the Hcrmi- 
tian quadratic form through the Heisenberg operator. 
Each Hermitian quadratic form has a dynamic matrix 
of its own. 

The Bogoliubov-Valatinian or Diracian diagonalizabil- 
ity of a Hermitian quadratic form is equivalent to the 
physical diagonalizability of its dynamic matrix. That 
is to say, the diagonalization of a Hermitian quadratic 
form is essentially an eigenvalue problem of its dynamic 
matrix. 

The dynamic matrix is always physically diagonaliz- 
able for a fermionic form. It may or may not be physi- 
cally diagonalizable for a bosonic form. Accordingly, the 
diagonalization exists and is unique for a fermionic form, 
forever. It exists and is unique for a bosonic form only if 
the dynamic matrix is physically diagonalizable. 

The dynamic matrix is the generator of the 
Bogoliubov-Valatin or Dirac transformation. The trans- 
formation required for diagonalization can be constructed 
immediately from the complete set of the orthonormal 
eigenvectors of the dynamic matrix, according to a stan- 
dard algebraic procedure. 

In a word, the eigenvalue problem of the dynamic 
matrix determines the diagonalizability of a Hermitian 
quadratic form, definitely and completely. 

Finally, it is worth emphasizing that the quadratic 
Hamiltonian is just regarded as a Hermitian quadratic 
form in this review, i.e., only its mathematical proper- 
ties are considered here. The physical instability as well 
as phase transitions of a system has not been concerned 
at all. That in itself is another intriguing problem, please 
refer to Ref. fsl. 
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